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I. Columbia University, Geography Department, Remote Sensing 

Project. 

f Since June 1975, Columbia University, Department of 
Geography has been conducting research on various aspects of 
applying Landsat and other remote sensor data to a broad 
range of topics related to agriculture, agricultural develop- 
ment, and the environment. The research is carried out in 

collaboration with the NASA/Goddard Institute for Space 

Studies (GISS) and is in support of the Institute's program 
in Earth Resources. 

The primary purpose of the research is to aid in develop- 
ing a global agricultural monitoring system which would rely 
on satellite data collection systems as major data sources. 
This involves the study of three specific topics: 1) agri- 

cultural change, 2) traditional agricultural systems and 
3) relationships between agriculture and the physical environ 
ment, i.e. the agroecosystem. 

Primary data sources for these studies are simulated 
advanced mul ti spectral sensors of various resolutions, 24- 
Channel mul ti spectral scanner (MSS) data, Landsat data, as 
well as other remote sensor data made available through the 
offices of GISS. MSS, Landsat, Skylab and aircraft spectro- 

l 

meter data are analyzed by means of GISS-developed software 
and computing facilities. 


1 


The faculty of the Department of Geography, Columbia University 
works closely with GISS personnel as well as supervises the student 
research assistants who may undertake individual projects within 
the framework of this proposal. 

Participating Columbia University faculty members from the De- 
partment of Geography are: 

Professor Kempton E. Webb, Chairman 
Professor Leonard Zobler 
Professor Colin J. High 
Dr. Jerry C. Coiner 

In addition to the four principal investigators, the staff 
includes five funded graduate students, three part-time student 
U technicians and four graduate researchers funded from other sources. 

The following papers were prepared by the staff during the 
reporting period. 

1. Neary, Patricia. "Monitoring Agricultural Transition in 
Southern New Jersey Using Landsat Data." M . A . Thesis, Depart- 
ment of Geography, Columbia University, January, 1978. 

2. Neary, Patricia J. and Jerry C. Coiner. "Monitoring Agri- 
cultural Systems," (to be presented as a Poster Session at the 
12th International Symposium on Remote Sensing of the Environ- 
ment, Manila, April, 1978). 

3. Coiner, J. C. and R. C. Bruce. "Use of Change Detection 
in Assessing Development Plans: A Philippine Example," (to 
be included in the Proceedings of the 12th International Sym- 

| posiLnn on Rem ote Sensing of the Environment, Manila, April 

■ 1978 ). 
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II. Experiments in Progress 


A. Support to Airborne Spectroradi ometer Studies 

GISS/Col urnbi a has developed a high resolution spectro- 
radiometer. This instrument acquires line trace spectra over 
a range from .43 to I micron for sequential areas approximately 
18 meters square. To aid in interpreting these spectra and 
relating them to specific ground locations a 35 mm photograph 
is acquired with every tenth spectrum.* 

Analysis of data acquired by this system over Imperial 

Valley, California has pointed to the existence of a spectral 

shift in the chlorophyll red absorption edge associated with the 

2 

heading of wheat and grain sorghum. 

Preliminary results of a study being conducted with data 
acquired over the LACIE Intensive Test Site in Finney County, 
Kansas further substantiate this red-shift phenomenon. Figure 
1 shows the ratios (averaged over each individual agricultural 
field) of data from two bands, 780-790 nm and 740-750 nm, of 
the GISS/Col umbi a spectroradi ometer . Fields in the corn and 
grain sorghum categories were planted within the month preceding 
the date of data acquisition (May 21, 1975); therefore, these 
categories consist primarily of exposed soil surfaces. The 


*Chiu, Hong-Yee and William Collins, "A Spectroradi ometer for Airborne 
Remote Sensing." Photogrammetric Engineering and Remote Sensing , Vol. 44 (4), 
April 1978, pp. 507-517. 

^William Collins, "Remote Sensing of Crop Type and Maturity," Photo - _ 
qrammetric Engineering and Remote Sensing , Vol. 44 (1), Jan. 1978, pp. 43-55. 






• FIGURE 1: AIRBORNE ^ECTRORADKMER FIELD AVERAGE RATIOS 



alfalfa fields are mature with predominately full canopy cover. 
The winter wheat is also mature, having been planted in September 
and October the previous year. As Fig. 1 indicates, mature, 
irrigated wheat is discriminated from the other crops and field 
conditions by ratios greater than approximately 1.10. Some con- 
fusion occurs with 3 wheat fields below the 1.10 ratio, 1 corn 
field and 2 alfalfa fields above. These anomalies cannot be 
accounted for with the available ground truth information. 

To corroborate the spectroradi ometer data, JSC helicopter 
FSS (S-191 H) data, acquired over the Finney County Test Site 
on the same day (May 21, 1975), was also analyzed. The ratios 
of two 20-nm bands (780-800 nm, 740-760 nm) were averaged for 
each field flown within the site and plotted by crop type and 
maturity (Figure 2). As with the spectroradi ometer data, mature 
(heading and fully headed) wheat are separable from other crops 
and field conditions by a band ratio value of approximately 
1.124. One alfalfa and 3 wheat fields are misassigned. 

Data from both instruments confirm the existence of the 
widening absorption edge initially identified by Collins in the 
Imperial Valley. Figure 3 shows the correlation between data 
acquired by each instrument for fields that were jointly covered 
by the two missions. The r 2 value for the nine jointly covered 
fields is .985. These results support the findings that both 
instruments are detecting the same phenomenon. 

Additional analysis is continuing on FSS data acquired over 
this site on an earlier date (May 14, 1975). 
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FIGURE 2: SPECTRAL BAND RATIOS FOR VARYING CROP STAGES 

FSS - (780-800nm)/(71K)-750nm) 

Fliie COUNTY, KANSAS HAY 21, 1975 
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FIGURE 2: SPECTRAL BAND RATIOS FOR VARYING CROP STAGES 

FSS - (780-800nm)/(7AO-750nm) 

FINNEY COUNTY, KANSAS MAY 21, 1975 
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B. Landsat Applications to Agricultural and Land Use 
Change Studies 

1. Monitoring Agricultural Transition in Southern New Jersey 

f 

This research, part of a broader effort to explore the 
possibilities of monitoring agricultural change with Landsat data, 

is outlined below and described more fully in the Appendix to this 
report. 


1.1. Background 
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For this study a 13,675 hectare area in southern New Jersey 
was chosen, an area which is known to be under considerable pressure 
to adjust to new agricultural circumstances. The study area is 
located about 40 kilometers south of Philadelphia, Pennsylvania, 
and 15 kilometers west of Vineland, New Jersey. 

Agriculture and food processing combine to form New Jersey's 
largest industry. Over the past few years, the state's total area 
in vegetables has declined whereas total area in wheat, corn, and 
soybeans has increased. Table 1 shows a decline in hectares har- 
vested over the period 1971-1976 for commercial vegetables, in 
particular, for processing vegetables. 1 


Significant discrepancies have been noted with respect to data in this 
category as reported annually in New Jersey Agricultural Stati^Hrs for the 
years 1971 through 1976. (Each aHHUaTFip fe-f Ogives statist! cs~ THThis category 
for preceding years as well as the current year.) For example the 1972 Statistics 
reported the total area harvested in 1971 for all commerce! vegetables TbF pro - 
cessing as 17,200 hectares. The 1976 Statistics reported 17.600 hprtaroc -in P +h-|c 

Tr 19 a ; d ii! ‘S 77 V 16 fi8 « Siven was 13,400 hectares, as si iil 
Table 1. According to the New Jersey Crop Reporting Service, the latest (1977) 
figures are the most accurate, implying that it took six years to rectify a 4,000- 
hectare error--an error amounting to nearly one-third the acreage in that cafponrv 
It appears that LANDSAT acreage estimators v/ould be more efficient in this regard.' 
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TABLE 1 


AREA HARVESTED IN HECTARES— ALL COMMERCIAL VEGETABLES 

STATE OF NEW JERSEY 
1971-1976 - 



Year 

For Fresh 
Market 

For 

Processing 

Total 


1571 

24,140 

13,400 

37,540 

1972 

24,120 

12,300 

36,420 

1973 

22,610 

12,970 

35,530 

1974 

20,850 

13,800 

34,650 

1975 

19,280 

11,080 

30,360 

1976 

19,850 

5,700 

25,550 


SOURCE: New Jersey, Department of Agriculture, New Jersey Crop 

Reporting Service, New Jersey Agricultural Statistics , August 1977. 


TABLE 1 


AREA HARVESTED IN HECTARES— -ALL COMMERCIAL VEGETABLES 

STATE OF NEW JERSEY 
1971-1976 . 



Year 

For Fresh 
Market 

For 

Processing 

Total 


1571 

24,140 

13,400 

37,540 

1972 

24,120 

12,300 

36,420 

1973 

22,610 

12,970 

35,530 

1974 

20,850 

13,800 

34,650 

1975 

19,280 

11,080 

30,350 

1976 

19,850 

5,700 

25,550 


SOURCE: New Jersey, Department of Agriculture, New Jersey Crop 

Reporting Service, New Jersey Agricultural Statistics , August 1977. 


Economic pressures have forced several vegetable processing concerns 
in New Jersey to move their operations elsewhere. Many New Jersey 
farmers have thus had to find either new markets or new crops. 

* The approximately 14,000 hectares which comprise the test 
area for this study lie in and around the area formerly farmed by 
Seabrook Farms, Incorporated, which closed its fresh vegetable pro- 
cessing plant in Cumberland County in March 1976. This meant that 
1,485 hectares farmed by Seabrook and 8,100 hectares under contract 
with 150 farmers would be released from their current land uses. 

All 1,485 hectares farmed by Seabrook are included in the test area 
except for a small area surrounding the town of Seabrook itself. 

Also included are privately owned farms, most of which operated in- 
depently of Seabrook; these are included as a comparison to the 
Seabrook land which had been completely devoted to vegetables. 

1.2. Experimental Design 

The purpose of this study is to detect agricultural systems 
change by separating two agricultural systems, corporate vegetable 
farming and field cropping, in order that the transition from the 
former to the latter may be detected. This study traces land 
utilization over time in an area where two agricultural systems 
exist in close proximity, each imposing a different spatial pattern 
on the landscape. 

With respect to "land utilization in time," crop calendars 
vary for each of the two types of agricultural system in this area. 
For example, vegetables are planted mainly during the warmer months, 
whereas winter grains are planted in the Fall. The period of active 
photosynthesis for vegetables extends later into the Fall than it 


10 


does for field crops. Winter grains are seeded in October and do 
not make an impact on land cover until early November. 

In addition to the temporal component of land utilization, 
the spatial component must be considered. In the agricultural area 
being examined here, vegetation and soil are the dominant surface 
features. Actively photosynthesizing plants are characterized by 
a relatively large amount of near infrared reflectance (bands 3 and 
4 on the MSS) and a relatively small amount of red reflectance (band 
2 on the MSS). Soils exhibit spectral patterns distinct from vege- 
tation while also differing within soil categories. Generally, 
spectral response curves of soils are rather flat in the visible 
region with a decrease in reflectance in the near infrared. The 
reflectivity of soils in the visible is also dependent upon soil 
color--a lighter soil will be more reflective than a darker one. 

The coastal plain soils in the test area are rather bright and there 
fore strongly influence spectral response in the visible bands. 

In this study, the ratio of band 4 to band 2 is used as an indicator 
of vegetation density and growth stage. Plotting the ratio of band 
4 to band 2 over time for a specific ground area yields a curve for 
each of the four years of the study, similar to the diagram below. 



TIME 
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Although Landsat MSS measurements of reflected energy are made 
at approximately the same time of day, the sun elevation at a specific 
latitude at a given hour changes throughout the year. To account for 
this variation, a sun elevation correction factor was introduced which 
adjusts all measurements to a solar angle of 90°. In addition, the 
use of the band 4/band 2 ratio reduces solar angle differences as 
well as signal fluctuations and differences attributable to changing 
atmospheric conditions by using a relative measure, a ratio. 

The analysis employed digital Landsat data and relied upon, 
as support materials, Landsat imagery, aerial photography, United 
States Geological Survey maps, Seabrook Farms Corporation maps, and 
ground-collected information. 

Although Landsat data are available for passes made every 
eighteen days from July 1972 to January 1975 and every nine days 
from January 1975 through most of 1976, not all of it is useable, 
due to excessive haze or cloud cover or to technical problems in 
producing computer compatible tapes of the data. Out of 119 pos- 
sible coverages, 26 acquisitions (or 22 percent) were found to be 
useable in this study. All 26 scenes were preprocessed for geo- 
metric correction and rectification to a Universal Transverse Mer- 
cator grid. 

Imagery was used initially to choose suitable passes and to 
geographically identify the corresponding digital data. A Landsat 
image, or “scene" corresponds to a block of data over a geographic 
area of approximately 34,225 square kilometers, and each data tape 
contains information from an entire scene. From each scene twenty- 
four smaller areas of 595 hectares, or "quads", are used as the 


basic unit of analysis. These areas may encompass up to 40 agri- 
cultural fields. It is felt that this unit is suitable for the 
long-run objective of the study--to simulate global monitoring, 
tn addition, the unit conforms to technical constraints in the soft- 
ware used in support of the project. Each set of twenty-four areas 
is registered spatially from date to date. Therefore the study 
concerns the change in reflectance of twenty-four identical areas 
over the course of four years--1973 through 1976. 

An unsupervised classification routine developed at GISS was 
used to obtain a nominal spectral response for each quad. Normally, 
this routine is used to delimit clusters of pixels which have similar 
response curves. In this study, however, the object was to delimit 
one single cluster for each quad into which approximately 95 percent 
of all the pixels in the quad fall. This percentage was chosen to 
filter out highly reflective, exotic pixels which might influence 
the quad's spectral response disproportionately. The classification 
routine was applied to twenty-four quads from each of twenty-six 
passes to obtain 624 characteristic 4-band spectral responses. The 
reflectances for bands 4 and 2 were ratioed for each of the 624 no- 
minal spectral responses and then plotted against time for each of 
the 24 quads. 

1.3. Results 

Classification was performed on all of the 24 quads for each 
of the available 26 LANDSAT passes and the band 4/band 2 ratios cal- 
culated and plotted. Two sample plots are presented in Figures 4a 
and b for quads 2 and 23, which are representative of field crops 
and vegetables respectively, as determined from information collected 


in the field and maps from Seabrook Farms. Figure 5 is a super- 

'*• imposition of these two plots to facilitate comparison. 

ILT ; 

The shapes of the curves differ from year to year because 
of, lack of concurrence in satellite acquisition times for 
usable data and local atmospheric differences. The two quads 
can be compared to each other within each year, however. The 
curves diverge within the summer months when both field crops 
and vegetables are growing. Field crops tend to be in full 
canopy during July and August whereas vegetables tend to have 
less canopy cover with rows being separated by strips of soil. 
The large amount of exposed soil decreases the band 4/band 2 
ratio, and this is evident in the figure. Curve separation may 
also be related to meteorological events which have immediately 
preceded data acquisition. For example, the Landsat pass during 
the summer of 1975 occurred ten to fifteen days after a series 
of severe storms destroyed crops in the southern New Jersey re- 
gion. In Figure 5, this event is reflected in the similarity 
of the two plots of 1975 during the summer, when they would be 


expected to diverge. 

The application of a non-parametric Mann - Whitney "U" test 
to the band 4/band 2 ratios of the two groups of quads (vegetable 
and field crops) yielded a significant difference in the means 
of the two groups of ratios for nine of the twenty-six acquisition 
dates, with a 98 percent level of confidence. Seven of these 
nine passes were during July and August, thus emphasizing the 
significance of the divergence between the two agricultural systems 

in the summer months. 
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Figure 4a Quad 2 • (field crops) 
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LANDSAT 4/ 2 Ratio of Vegetable and Corn/Soy Systems 

So. New Jersey, 1973-1976 



Vegetable System 
Corn /Soy System 




To test the separability of the two agricultural vegetation 
systems from the natural vegetation component of the landscape, 
the nonparametri c Kruska - Wallis test was applied to the three 
sample groups. This test revealed that the three sample groups 

f 

belong to different populations, at a 99 and 95 percent confidence 
level, respectively, for the data from June 1973 and June 1976. 

As a preliminary test of an aggregated-uni t approach to 
monitoring agricultural systems, the results are encouraging. The 
Landsat derived radiometric indicator, even at this level of 
aggregation, appears to offer a viable method of differentiating 
agricultural systems from each other as well as from natural ve- 
getation elements. 

2. East African Swidden Agriculture 

An investigation into the feasibility of detecting shifting 
cultivation in East Africa has encouraging preliminary results. 

The study site is in the West Pokot District of Kenya. Landsat 
computer-compatible tapes for three dates are available: scene 
1067-07221, 9 September 1972; scene 1193-07230, 1 February 1973; 
and scene 2063-07111, 26 March 1975. 

Supplementary data include 1956 black-and-white aerial photo- 
graphy and data from 1961-62 anthropological fieldwork. 

Using the supplementary data in conjunction with gray scale 
printouts of geometrically corrected LANDSAT data, a preliminary 
study site of 595 hectares was selected. An approximately ten- 
hectare area in this study site lies in a region locally known 


as Asar and always contains active cultivation. This area was 
thus selected as the training area, from which the most frequent, 
or modal, response in each spectral band was determined, for each 
of the three available dates. These modal responses constituted 
the training data. 

These training data, after successful application to the 
preliminary study site, were extended to a larger 200 by 200 
pixel area (approximately 19,000 hectares). 

Information about the study area and swidden agriculture as 
it is practiced there, in conjunction with an understanding of 
the characteristics of the satellite data and of the analysis 
procedure, provides a basis for the development of a model of 
expected results. The model developed here has three components: 
data points classified as swidden should be on east-facing slopes 
below 2,500 meters; one-sixth to one-fourth of the cultivable 
east slope area below 2,500 meters could be in cultivation at 
one time; a number of data points should be in swidden on all 
three dates, covering two and one-half years. Analysis results 
were evaluated in terms of these three expectations: for the 
three dates, 78.6 percent, 88.9 percent, and 86.5 percent of 
the points classified did fall on east-facing slopes below 2,500 
meters; roughly one-eight of the cultivable east slope area below 
2,500 meters was in cultivation on each date; 387 data points 
were in swidden on the first two dates, and 76 data points were 
in swidden on all three dates. 

The high rate of correspondence between expected results and 
analysis results leads to the inference that swidden agriculture 


has been successfully identified in this area using the method 
described above. 

3. Suburban Land Use Change and Flooding 

Major flooding has occurred in the Saw Mill River watershed, 
Westchester County, New York, in 1975 and 1977 despite the con- 
struction of numerous flood control structures. It is believed 
that watershed development is a major contributing factor to the 
flooding, due to reductions in natural vegetation, compaction of 
the earth's surface, increases in impervious surface area and in 
smoother less frictional surfaces. The purpose of this study is 
to monitor changes in land use categories of hydrologic signifi- 
cance over a five year period (1972-1977) using Landsat data. 

The areal extent of the following land use categories, each 
of which posesses unique runoff characteristics, will be deter- 
mined for the 26-square mile watershed for both 1972 and 1977; 

1) Commercial, business and industrial 

2) Residential- multi-family housing 

3) Residential- single family homes situated on % acre lots 

or smaller 

4) Residential- single family homes situated on lots greater 

than H acre 

5) Lawns, parks and cemetaries 

6) Woodlands 

7) Water 

Aerial photography (1:24,000) acquired in 1970 and 1976, in 
addition to ground survey of the area, will be used to assist in 
the selection of test sites for the development of spectral sig- 
natures for the land use categories. 


C. Effect of Environmental and Agro-Management Factors on 
Computer Classification of Crop Type 

Current research in computer classification of remote sensing 
techniques concentrates on acreage estimation of agricultural 
land use. The techniques used for acreage estimation identify 
and then separate radiance values of the scene into categories. 
These values are usually considered discrete elements of the en- 
vironment, and no attention is given to their composition. Another 
way of examining the spectral response of any landscape is to 
view it as a composite of management and physical environmental 
characteristics, which make varying contributions to the land- 
scape's radiance. It is the purpose of this research to assemble 
major environmental and agro-management factors into a model 
measuring their relative influence on spectral signatures of an 
agricultural landscape. The area of study is in a LAC I E inten- 
sive study site in the northern Great Plains. The remote sensor 
data was acquired by the Johnson Space Center 24-Channel multi- 
spectral scanner, August 15, 1975. 

The computer classification method used in this digital 
analysis consists of two major elements: 1) the ground truth 

system and 2) the radiance values of the landscape. The first 
element divides the scene into finite geographic units of varying 
sizes and shapes (i.e., agricultural fields, soil boundaries), 
and the second element represents the spectral response of the 
environment. It is the function of the classification algorithm 
to distinguish between and assign each pixel to its respective 
land use/land cover category. A change in the ground truth system 
results in classification of the same landscape into different 


categories. This may reveal a greater degree of homogeneity 
within a category, resulting in higher classification accuracies 
and a more thorough understanding of the spectral reflectances. 

In agricultural studies, two categories of environmental 
parameters are Identified: management and physical environment. 
Management consists of such farming factors as tillage direction, 
fertilization, planting date, and irrigation. The physical en- 
vironment category includes elements such as soil type, soil color, 
and terrain slope. The relative influence of a set of parameters 
on land use classification is unique for each landscape. 

The study region for this investigation is a 5 x 7.S kilometer 
area of Williams County, North Dakota, where severe winters and 
dry summers predominate. Spring wheat is the major crop in a 
landscape intermixed witii'summer fallow and pasture fields. 

Summer fallow is practiced almost exclusively as dn anti-wind 
erosion measure in strip cropping patterns. Pasture includes 
harvested and non-harvested oats and other natural grasses which 
may or may not be grazed. The agricultural practices are in- 
fluenced to some extent by remnant glacial features, such as 
potholes and till knobs. Potholes create areas of poor drainage 
which become marshy during stages of evaporation after the rainy 
season. Sections of patchy bare soil (till knobs) produce a 
mottled effect in the appearance of the landscape. 

The management practices and physical environmental factors 
of this area comprise a landscape system of varying spectral res- 
ponses. Explanation of the spectral variability will be attempted 
through the testing initially of three environmental /management 
parameters as separate ground truth systems: tillage direction. 


planting date, and soil type. 


To date, the three ground truth files have been constructed 
In preparation for computer classification of the mul ti -spectral 
scanner data. The ground truth management system used in this 
study is a computerized geographic information system developed 
in cooperation with GISS.* 


Continuing work on this experiment involves computer classi- 
fication using each of the three ground truth systems in order 
to determine the relative influence of each environmental and 
agro-management parameter on classification accuracy for the 
land cover classification categories. The classification pro- 
cedure is as follows: 

1) Flight Line 1 was selected as the training area. An 
unsupervised classification is performed on this data in 
order to derive spectral signatures for each individual 

' land cover category. The ground truth arrays are used to 
"mask" the data. For example, to determine the signature 
of the "wheat-soil type A" category, only pixels in that 
category, as defined by the ground truth, are used in the 
classification routine. 

2) Once signatures are established for each combination 
of crop type and management or envionmental parameter 
value, all signatures for land cover categories in a single 
management or environmental parameter ( e . g . , soil type) are 
simultaneously applied to the entire training area in a su- 
pervised classification mode to determine classification 
decision criteria. Flight lines 2 and 3 are then classified 
using these signatures and their associated decision criteria 
as established from the training data. 


■^See Coiner, J. C. and Ungar, S. G., "Ground Truth Management System to 
Support Multispectral Scanner (MSS) Digital Analysis." Proceedings of the 
ASP-ASCM Joint Annual Meeting, Washington, D. C., February 27 - March 5, 1977, 
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INTRODUCTION 


Purpose 

Agriculture is one of this nation's most important 
industries. It is also one of the prime activities through 
which man modifies the natural ecosystems of the earth. It is 
therefore of tremendous importance that decision makers in 
business and government have an accurate picture of the state 
of agricultural production from season to season and of the 
long-range trends in different agricultural regions. This 
paper will explore the possibilities of monitoring agricul- 
tural systems change with LANDSAT data and will discuss the 
results of one attempt using data from southern New Jersey. 

Why Monitor Regional Agriculture? 

The ability to detect and measure agricultural change 
allows a geographer to assess the impact of new land uses 
on the social and business community and on the environment 
of a particular region. There are a number of reasons for 
investigating the type and amount of agricultural change 
which may be. occurring. First, there are the economic con- 
considerations. Local agribusiness must identify regional 
trends in order to anticipate future demands of specialized 
services. A fertilizer manufacturer must know what crops 
are gaining in importance so that he can plan the production, 
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distribution, and storage of his product. A dealer in farm 
machinery would be interested in agricultural trends to aid 
in planning for future machine inventory and maintenance 
demands. Builders of storage and marketing facilities 
monitor regions undergoing conversion of cropping systems to 
enable them to provide necessary warehousing. Commodity 
brokers, wholesalers , and retailers have an interest in how 
agricultural change will affect prices. Local governments 
use agricultural statistics to establish assessment values 

and tax rates. 

Secondly, environmental conditions are affected by 
agricultural change. Different cropping systems require 
different types of fertilizer, herbicide, and pesticide which 
affect the environment when used and, perhaps, for a long 
time afterward. Irrigation practices can affect groundwater 
supplies or alter flow in streams. In addition, certain crops 
can encourage or discourage organisms from entering a parti- 
cular environment. 


1 Hill-Rowley et al., in -“Improved Resource Use Decisions 
and Actions through Remote Sensing," National A eronautics an d 
Space Administration Earth Resources Survey Sy mposium (Houston • 
National Aeronautics and Space Administration, June ±3/b), 
1747-1768, note that for economic and technical effective- 
ness processing facilities for beans, corn and small grains 
should be located very close to areas of high production. 
Specifically their report states ". . .in the past, elevator 

location decisions were a 'hit-or-miss' process based main y 
on subjective and unscientific procedures and criteria, crop 
Reporting Service statistics regularly indicate production 
totals by counties. However, elevator service areas, with a 
radius of 13 to 16 kilometers, are areas substantially smaller 
than a county, and many include parts of more than one county. 
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LANDSAT as a Monitoring Tool 
Since agriculture represents a renewable resource 
whose cycle of renewability can be as short as two months, 
monitoring the resource is an onerous and expensive task. 
Traditional methods of data collection such as ground 
surveys and aerial photographs are more appropriate for 
slowly changing phenomena, for example, nonrenewable resources, 
geologic features, or certain residential areas. "Monitoring" 
has been described as the process of "obtaining a new inven- 
tory periodically .” 1 LANDSATs I and XI, earth resources 
satellites maintained by the National Aeronautics and Space 
Administration, offer a means to collect data periodically 
at little expense to the user. Between 1972 and 1975, a 
LANDSAT satellite passed over the same locale every eighteen 
days and from January 1975 through most of 1976, every nine 
days. It is, therefore, well suited for monitoring a seasonal 

resource such as agriculture. 

Were LANDSAT used to monitor regional agriculture, it 
could prove very economical for the agribusiness community. 

A recent study estimates that the improved information capa- 
bilities afforded by LANDSAT could benefit the United States 
annually a sum of $200 to $250 million for wheat, $50 to 


'•Robert N. Colwell, "introduction,” in Mon|| l of Rem ote 
Sensinq, ed. Robert G. Reeves (Falls Churchy Virginia. 
AmerlS m Society of Photogrammetry , 19 75), p. 2. 
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S1 00 million for corn, and $6 to $11 million for soybeans, 
"using conservative assumptions on expected LANDSAT system 

..1 

performance. 

The methodology proposed in this paper is tested 
area in southern New Jersey, where it is known that a change 
in land use is occurring due to the loss of a specific 
market. An indicator is developed which separates two 
agricultural systems in order that the change from one to 

t. ^ should the monitoring of agricul- 

the other may be detected, bhouxo t 

. . • mau iprsev prove feasible/ it may be 

tural transition m New Jersey piuvu 

■ nunr.i reci ions which arc currently 
extended to other agricultural regions 

undergoing rapid change. 


1 Econ, incorporated, EcowmicBeng | l t8 
^^^^fXc^HnH^ScHnX^alTn^^ma^Ion N77-22003, 

NASA-CR-152657, 1977), p. ilx * 



AREA OF INVESTIGATION 


Location and General Characteristics 
For this experiment a 13,675 hectare (approximately 
34,000 acres) area in southern New Jersey was chosen, an 
area which is known to be under considerable pressure to 
adjust to new agricultural circumstances. Map 1 shows the 
approximate location of the study area, which is located about 
40 kilometers south of Philadelphia, Pennsylvania, and 15 
kilometers west of Vineland, New Jersey. Map 2 shows more 
specifically the area which is being analyzed. 

This specific region lies in both Salem and Cumber- 
land counties, within the Upper Pittsgrove and Alloway 
townships in Salem County and the Upper Deerfield Township 
in Cumberland County. Upper Pittsgrove and Alloway are in the 
central portion of Salem County. The landscape is gently 
rolling, underlain by geologic formations of loosely con- 
solidated sedimentary deposits. Soils in the parts of these 
townships which are in our test area belong to a soil asso- 
ciation that is silty, well-drained to poorly drained, nearly 
level to gently sloping. 1 In 1953 this part of Salem County 

1 U.S., Department of Agriculture, Soil Conservation 
Service, New Jersey Agricultural Experiment Station, Soil 
Survey of Salem County, New Jersey , by Van R. Powley 
(Washington , D.C. : Government Printing Office, 1969) p. 5. 


of Test Site within Cumberland 


ORIGINAL PAGE IS 
OF POOR QUALITY 


SOURCE: United States Geolcgicao. Survey 

1:250,000 map, entitled " Wilmington , Delaware, 
New Jersey; Penna. ; Maryland." 
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was characterized by "dairy-cash crop and corporation vegetable 
farming on heavy textured soils. Personal observations in 
1976 reveal that the Salem County section of the test site is 
dominated by corn, wheat, soybeans, potatoes, orchards, 
pastures, and nurseries. 

Upper Deerfield Township in Cumberland County is immedi- 
ately south of the area just described. Its topography and 
soils share the same general properties. Most of that part 
of the township which is in the study area is very well suited 

for agriculture. In 1974, 5,158 hectares (12,746 acres) in 

2 

the Upper Deerfield were devoted to agriculture. Until 1976 
vegetables dominated the agriculture in the township; beginning in 
1976 corn and soybeans replaced vegetables as the dominant crop. 

Recent Agricultural History 

Agriculture and food processing combine to form New 
Jersey's largest industry. The position of the state in the 
Boston-to-V7ashington urban sprawl along the east coast of the 
United States affects its agriculture both positively and 
negatively. The advantages of its position are obvious in 
that it lies in close proximity to urban centers which serve 

"^Leonard Zobler, "Man-Land Relations in Salem County, 

New Jersey" (Ph.D. dissertation, Columbia University, 1953) , 
p. 193. 

^Cumberland County Planning Board, Land Use: Existing 

Patterns and Environmental Characteristics, Upper Deerfield 
Township. Cumberland Co unty, New Jersey (Bridgeton, New Jersey; 
Upper Deerfield Township Planning Board, 1974), p. 35. 

^Donald Janson, "State to Gauge Pollution in Hope of 
Easing Curbs," New York Times , 16 January 1976, p. 63. 
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as its markets for both fresh and processed foods. The state 
has been known for many years for its "commercial, gardenrng 
and fruit culture." 1 On the other hand, the nearness to 
urban centers increases the problems of taxation, labor, 
development, and environment, many of which result from compe- 

. v-o for a comparable standard of living, 
tition with urban centers for a comparau 

. . u- wo effects on New Jersey agriculture: 

These problems have had two etrecrs 

, - . _ _ n a t he sellout to industry and hous- 

(1) the abandonment of farming and the seuo 

, /0 i _ aradual change from labor-intensive 

ing development, and (2) a gradual cn j 

,, -Fvri-iiO to more mechanized 

systems (for example, vegetables and fruxts) 

farming (such as wheat, corn, and soybeans). 

Over the past few years, the state's total area in vege- 
tables has declined whereas total area in wheat, corn, and 
soybeans, has increased. 2 The area harvested in Cumberland 
and Salem counties for selected field crops is detaxled xn 
Table 1. These latter crops have always played a part xn the 
New Jersey agricultural system, but their recent areal xncrease 
is indicative of a statewide trend. Table 2 shows the declxne 


1 This type of farm! ng is described as ” h J?^ s “Hoir 

^fs^ n uy fe^Sd 1 ?^ ^“—r^ srs&r 

do not require lar |® “ "““^““ith ed. (Hew York: McGraw Hxll 
et al.. Elements of Geography, 

Book Company, 1957), p. 591. 

2 New Jersey, Department of Agr icuttui e^NewJer seyPf 
Reporting sSvic4, 

1977. 
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TABLE 1 


HECTARES — SELECTED FIELD CROPS 
AREA harvested in hectare x 

CUMBERLAND AND SALEM COUNTIES, 1971-1976 



Cumberland 

1971 

1972 

1973 

1974 

1975 

1976 


Salem 

1971 

1972 

1973 

1974 

1975 

1976 


1,290 

1,150 

1,320 

1,340 

1,460 

2,020 


930 

1,290 

1,620 

1,940 

2,270 

2,350 


5,260 

4,130 

4,740 

5,020 

5,220 

6,480 


1,340 

1,660 

1,990 

3,360 

3,160 

3,320 


1,170 

1,170 

1,340 

1,700 

1,900 

3,040 


2,020 

2,430 

2,910 

3,640 

3,850 

5,950 


SOOTCE: New 


2, August. 

Vne area harvested in hectares f has ^been^calculated h * ctares 

author, using the ^ n ^^°^rves?ed provided in the NJDA sources, 
from information on acres n 


TABLE 2 


AREA HARVESTED IN HECTARES — ALL COMMERCIAL VEGETABLES 

STATE OF NEW JERSEY 
1971-1976 1 


Year 

For Fresh 
Market 

For 

Processing 

Total 2 

1971 

24,140 

13,400 

37,540 

1972 

24,120 

12,300 

36,420 

1973 

22,610 

12,970 

35,580 

1974 

20,850 

13,800 

34,650 

1975 

19,280 

11,080 

30,360 

1976 

19,850 

5,700 

25,550 


SOURCE: New Jersey, Department of Agriculture, New Jersey Crop 

Reporting Service, New Jersey Agricultural Statistics , August 1977. 

^The area harvested in hectares has been calculated by the author, 
using the conversion factor of 1 acre = .40468564 hectares, from infor- 
mation on acres harvested proved in the NJDA sources. 

2 The amount listed in the "Total" column is the sum of hectares 
listed in the previous two columns. 
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in hectares harvested for commercial vegetables, both fresh 
market and processing.^ Prom this table it appears that the 
processing vegetables have undergone the more severe drop. 
Economic and environmental pressures have forced many members 
of the vegetable processing community in New Jersey to move 
their operations elsewhere. In recent years, the Del Monte Cor- 
poration, P. F. Ritter (Curtis-Burns Corporation), and the 
Seabrook Foods Corporation have ceased purchasing New Jersey- 
grown produce. 2 New Jersey farmers have had either to find 
new markets or new crops; if they choose the latter alternative, 
field crops are often the most suitable because (1) a farmer 
can easily convert to this labor-extensive system and (2) wheat, 
soybeans, and corn have been drawing an increasingly higher 
price on the market (wheat and corn dropped in 1977) . 

The approximately 14,000 hectares which comprise the 
test area for this study lie in and around the area formerly 
farmed by Seabrook Farms, Incorporated, which closed its 


^These figures were also obtained from New Jersey . 
Agricultur al Statistics (1977) . It was noted that there was 
a significant discrepancy between data reported for the years 
1971 through 1976 in each year's report. For example New 
Jersey Ag ricultural Statistics (1972) reports that acres 
hirFvested for all commericial vegetables for processing amounted 
to 42,420. In 1976 Statistics reports 43,420, and in 19// it 
reports 33,120. That the latest figures are supposed to be 
the roost accurate (according to the New Jersey Crop Reporting 
Service) implies that it took six years to find a 10, 000-acre 
error, an error amounting to one-third the acreage m that 
category. It seems that 1ANDSAT acreage estimators would be 
more efficient in this regard. 


2 Donald Jansen, "Seabrook Farms' Plan to End Processing 
of Vegetables is Termed a 'Disaster,'" New York Times, 

6 January 1976, p. 67. 
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fresh vegetable processing plant in Cumberland County in March 
1976.^ This meant that 1,485 hectares (3,670 acres) Of their 
own and 8,100 hectares (20,000 acres) under contract with 150 
farmers would be released from their current land uses. All 
1,485 hectares farmed by Seabrook are included in the test area 
except for a small area surrounding the town of Seabrook itself. 
Also included are privately owned farms, most of which operated 
independently of Seabrook; these are included as a comparison 
to the Seabrook land which had been completely devoted to 
vegetables. 

The departure of one processor does not ordinarily sig- 
nify a trend; however, in this case the reasons given for 
such action relate to regional and state problems and are 
similar in nature to those expressed by Ritter and Del Monte. 
The president of Seabrook Farms, James M. Seabrook, cited the 
following causes for the January 1976 decision: 

a) The fringe benefits to labor in New Jersey are higher 
than in the other states where the company operates 

b) The farmers' prices for the growth and delivery 
of vegetables are higher in New Jersey than elsewhere 

c) The environmental requirements posed by state law are 
difficult to implement; for example, the Department of Environ- 
mental Protection refused a request to burn high-sulphur fuel 
oil, a measure which Seabrook and others considered necessary 

repackaging plant and warehouse are still maintained 
by the Seabrook Farms Division of Seabrook Foods, Atlanta, 
Georgia, a wholly-owned subsidiary of Springs Mills, Inc. The 
plant and warehouse now handle only bulk frozen produce, 
however, which has been grown and packed elsewhere. 


economically 


d) New Jersey is less pioneering in its research into 

vegetable problems than California and other states 1 

2 

e) energy is less costly in other states. 

It appears, then, that the Seabrook action, along with 
those of Ritter and Del Monte before it, are part of a larger 
agricultural change occurring in the state. As James M. 
Seabrook commented on the situation in January 1976, "The 

T « 3 

future is bleak for vegetable farmers m New Jersey. 


Ws point is debatable. C. W. Thornthwaite did pio- 
neering work in microclimatology and irrigation ° f e^^for 
in conjunction with Seabrook Farms m the 1950s. Si T°^ tion 
example, C. W. Thornthwaite, Climate and Scientific Irrigatio n 
in New Jersey, The Johns Hopkins University Laboratory of 
Climatology Publications in Publications in Cl ^™ atol °^ 
vol. 6, no. 1 (Seabrook, New Jersey, 1953 ) . Also, see wor 
hv others both in the 1950s and more recently, such as J. J. 
nLni ns/ instructions for Making Phenological Observations 
of Garden’T^T'The Johns Hopkins university Lab ° r ^^^ v Cl1 a 
tology Publications in Climatology, vol 5, ^ (S^brook, 

New Jersey, 1952); David K. Robertson, Gr ound Water Ava ^La 
bi litv in Southern New Jersey , C. W. Thornthwaite Associates 

Publications in Climatology, vol. 26, no. 1 N ®" 

Tersev 1973) ; and J. R. Mather, Estimation of Areal Average 
S i “”:'.i s ka Different Network Dens ities and AverSjlSg. 
Techniques, C. W. Thornthwaite Associates Publications m 
ClimatoTogy , vol. 28, no. 2 (Elmer, New Jersey, 1975). 

2 Janson, "Disaster," p. 67. 


Ibid. 
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This being the case. it is imperative that monitoring of 
Hew Jersey agriculture he undertaken so that the sooiai ana 
business communities are prepared to accommodate the transr- 

tion . Traditional methods must he relied upon for now, hut 

1 out the potential usefulness of 

the New Jersey example porn 

Modern data collection tools, 

a fast and inexpensive method. Modern e. 

such as the Multispectral Scanner described below, are use 
in this study to develop such as the Multispectral Scanner 
described below, are used in this study to develop such a 


method . 
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DESCRIPTION OF THE LAND SAT SENSOR AND 
THE DATA IT OFFERS 

LANDSAT I and II serve as platforms for a particular 
remote sensing device , 1 the multispectral scanner (MSS). 

The MSS consists of four sets of detectors which scan the ^ 
earth as the satellite passes in its sun-synchronous orbit. 

Each set of detectors is sensitive to a different portion of 
the solar radiation reflected from the ground object. 

Table 3 lists the four bands of the MSS and the portions of 
the electromagnetic spectrum to which each is sensitive. 
Recipients of energy from the sun may reflect, transmit, or 
absorb (and then re-emit) the sun's energy. This applies to 
both targets on the earth's surface and particles in the 
atmosphere. The MSS measures mainly the reflected energy 
from both these recipients . 3 

1 0ther instruments use LANDSAT as a platform but are 
not discussed here. 

2„ q synchronous" is defined as "as an Earth satellite 
orbit in which the orbital plane ^““^“^"placls on 
altitude such that the sa _ twice daily at the same local 
s^UmlT-^n^meScan Society of Photogrammetry , 

Remote Sensing , pp. 2105-2106. 

Wing daylight hours, gSfigg “^fnricro- 

radiation returned to space begins to dominate, according 

meters, w h f l ^S'^^ i ^^°5 h otographic and Nonphotographic 
to John Estes, Imaging wxc _ J John E. Estes and 

sensor Systems = ^"l^lfjIgl^^Son^lIEing Co., 
Leslie W. Senger (Santa Baioara. noiui 

1974) , p. 36. 


TABLE 3 


THE FOUR BANDS OF THE MSS , THEIR CORRESPONDING 
WAVELENGTH INTERVALS IN MICROMETERS AND 
THEIR COMMON NAMES 


MSS 

Wavelength 

Common 


Interval 

Name 

1 

0.5 - 0.6 

Green 

2 

0.6 - 0.7 

Red 

3 

0.7 - 0.8 

First Near- 
Infrared 

4 

0.8 - 1.1 

Second Near- 
Infrared 


LANDSAT data products arc of two general types, namely, 
imagery and digital tapes. The first, imagery, is a remote 
sensing term which includes both photographic and nonphoto- 
graphic (but photograph-like), products of sensing instruments, 
from cameras to radar. In this study we use imagery only as 
a locating tool. LANDSAT imagery is a pictorial display much 
like an aerial photograph but produced from digital data 
acquired one scan line at a time. Data represent amount of 
radiation reflected from the ground, and intervals of reflectance 
are assigned gray tones depending on this amount. In this 
study, imagery is used as a device to associate portions of 
the digital data with known ground areas. Digital information 
displayed as imagery allows one to use the shape and context 
of certain parts of the data to ascertain the ground features 
with which the data are associated. For example, if one 
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intends to analyze a ground area just east of a river, one 
would first examine imagery to determine the location of the 
river, easily discernible by the linear, dendritic pattern 
formed by the riverine vegetation which borders it. One would 
then know the approximate coordinates of the area on the 
digital tape. Digital tapes are our main source of data 
for analysis. 

The word "digital" here implies numbers assigned to 
physical entities. In LANDSAT data the physical entity with 
which one works is a parcel of ground surface. The smallest 
ground area for which reflection is measured — that ground area 
to which a number is assigned — is approximately .464 hectares 
or 60 by 76 meters on a side. This unit of ground area is 
called a "pixel" ("picture element") when referring to the 
data. The same pixel is represented in four bands; in other 
words, four measurements are made for each pixel, or .464 
hectares. See Figure 1 for a diagrammatic representation of 
this . 

The spectral response of a pixel is shown in bar graph 
form in Figure 2. This is important because the shape of 
a pixel's spectral response curve may be associated with a 
particular type of land cover. In addition, aggregations of 
pixels may be spatially related to the same land cover class 
(for example, the numerous pixels in a large agricultural 
field) and thus have similar spectral response curves. 
Generally, digital analysis of LANDSAT data in agriculture 
involves clustering pixels into areal units associated with 
different crop types. In this experiment, however, the 
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Fig. 1. Five adjacent pixels <jna 
responses in units of energy (mw/cra / s 
of the four bands of the MSs. 
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areal unit is not one associated with a homogeneous crop class, 
derived from pixel similarities, but a unit of larger dimensions 
composed of 1280 pixels. Such an aggregation may include many 
agricultural fields and several crop classes. The purpose 
of aggregating a number of fields into one unit is to facili- 
tate a large area monitoring operation— as for the globe— in 
which field-by-field analysis would be too time-consuming and 
expensive. The units used, here called "quads" (a software 
term used at the Institute for Space Studies, Earth Resources 
Program, to describe a 40 by 32 array of LANDSAT data) , consist 
of 1,280 pixels or 595 hectares. Figure 3 shows four gray 
scales (corresponding to four bands) for one quad of 1,280 
pixels in the test area. Gray scales, like imagery, are 
simple displays of digital data in which gray tones, achieved 
by character overprinting, are associated with discrete inter- 
vals of energy responses bright with high responses, dark 
with low responses. For example, the printed character in the 
upper left hand corner of each gray scale represent the 
radiance/response for the same ground parcel in each band. 

The gray scales are shown here to give a visual impression 
of the data and the quad size. 
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EXPERIMENTAL DESIGN 
Theory 

The purpose of the work described herein is to detect 
agricultural systems change by separating two agricultural 
systems, corporate vegetable farming and field cropping, in 
order that the transition from the former to the latter may 
be detected. DeSchlippe describes "agricultural systems" 
in terms of four characteristics, two of which serve as the 
basis for this study, namely: 

(1) land utilization in space (pattern of field 
types on their respective ecological back- 
ground) 

(2) land utilization in time (pseudorotations) 1 
(The other two characteristics involve labor and nutrition, 
which are not appropriate in a remote sensing context.) 

This experiment traces land utilization over time in area 
where two agricultural systems exist in close proximity, 
each imposing a different pattern on the soil background. 

With respect to "land utilization in time," crop 
calendars vary for each type of agricultural system. For 
example, vegetables are planted mainly during the warmer 
months, whereas winter grains are planted in the Fall. 

1 Pierre DeSchlippe, Shifting Cultivation in Afri ca 
(London: Routledge and Kegan Paul, 1956), p. 238. 
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Figures 4 and 5 show the planting and harvesting dates for 
both field crops and vegetables respectively. Table 4 shows 
the approximate periods of peak growth for a variety of crops, 
calculated from Figure 4 for field crops and from a list of 
planting and harvesting dates used at Seabrook Farms. It 
appears from this list that the actively photosynthesizing 
period for vegetables extends later into the Fall than it 
does for field crops. Winter grains are seeded in October 
and do not make an impact on land cover until early November. 
Vegetables and field crops have somewhat different cyclical 
patterns, and this method attempts to trace these cycles over 

time. 

With respect to "land utilization in space," field crops 
and vegetables impose different percentages of canopy cover 
on the land. The influence of exposed soil on a spectral 

response is discussed below. 

Surface features influence the amount of solar energy 
reflected into space from the earth to an earth-oriented sensor. 
Vegetation, soil, exposed rock, and water reflect energy 
differently in each spectral band. In the agricultural area 
under examination in this paper, vegetation and soil are 
the most dominant surface features. Although vegetation types 
differ in their spectral response curves, they do exhibit 
general characteristics distinct from other types of land cover 

■^Information received from Mr. Val Perri, Seabrook Farms, 
Inc., August 4, 1977. 





Figure 4. New Jersey: Field Crops, Usual Planting 

and Harvesting Dates. 

SOURCE: New Jersey, Department of Agr i cu 1 tu re , New 

Jersey Crop Reporting Service, New Jersey Agricu ltural, 
Statistics , August 1976. 
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Figure 5. New Jersey: Vegetables, Usual Planting 

and Harvesting Dates. 

SOURCE: New Jersey, Department of Agriculture, New 

Jersey Crop Reporting Service, New Jersey Agricultura l 
Statistics, August 1976. 
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TABLE 4. PEAK GROWTH PERIODS OF 
FIELD CROPS AND VEGETABLES 


Mar 


Apr 


Cauliflov7er 
Carrots 
Broccoli 
Spinach — Fall 

Spring 

Squash 
Kale — Fall 

Spring 
Baby Liiaas 
Green Beans 


X 

X 


X 


X 


May 


X 

X 

X 

X 


X 


Jun 


X 

X 


X 


X 


Jul 


X 

X 

X 


X 

X 



Actively photosynthesizing plants are characterized by a 
relatively large amount of near infrared reflectance (Bands 3 
and 4 on the MSS) and a relatively small amount of red reflec- 
tance (Band 2 on the MSS) . Soils exhibit spectral patterns 
distinct from vegetation although they also differ among them- 
selves. Generally spectral response curves of soils are rather 
flat in the visible region with a decrease in reflectance in 
the near infrared. The reflectivity of soils in the visible 
is also dependent upon soil color — a lighter soil will be more 
reflective than a darker one. The coastal plain soils in the 
test area are rather bright and will therefore be very influ- 
ential on the visible bands' spectral responses. 

The increase in near infrared reflectance which occurs 
as green plants develop is associated with the maturation of 
the plant leaf. Simply described, the structure of a plant 
leaf, from upper to lower layer, consists of cuticle, upper 
epidermis, mesophyll (palisade parenchyma, then spongy paren- 
chyma), lower epidermis, and cuticle. The maturation of plant 
leaves is accompanied by an increase in the number of inter- 
cellular air spaces in the spongy mesophyll layer. The cell 

^S. J. Kristof and M. F. Baumgardner, "Change of Multi- 
spectral Soil Patterns with Increasing Crop Canopy," LARS 
Information Note 102372, (West Lafayette, Indiana: The 

Laboratory for Applications of Remote Sensing, Purdue Univer- 
sity, 1972) , p. 16. 
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walls of the spongy parenchyma act as multiple reflectors of 

near infrared radiation. As Gates describes it, 

A small amount of light is reflected from the 
leaf cuticle; much is transmitted into the 
spongy mesophyll, where the rays have frequent 
encounters with the cell walls ... The multiple 
reflection is essentially a "random walk" as 
the rays frequently change direction within 
the leaf. 1 

The reduced amount of reflected radiation in the red 
band is due to the absorption of light in the .6 to .7 micro- 
meter region by chlorophyll-a, the major pigment involved in 
photosynthesis. Five photoreceptive structures are found in 
plants as follows: grana-containing chloroplasts, grana-free 

chloroplasts , chromatoplasm, chroma tophores , and free grana. 2 
The grana— containing chloroplasts are most common in higher 
plants and are commonly found in the palisade layer. (In 
plants where the mesophyll is undifferentiated, such as in 
corn and other grasses, cells containing chloroplasts are 
found in herbaceous stems, twigs, and unripe fruits. ) Light 
enters the leaf and certain portions of its spectrum are 


■^David M. Gates, "Physical and Physiological Properties 
of Plants," in Remote Sensing (Washington, D.C. : National 

Academy, 1970), pp. 226-227. 

2 J. B. Thomas, Primary Photoprocesses in Biology , 
(Amsterdam: North Holland Publishing Co. , 1965) , pp! 127-135. 

2 Carl L. Wilson and Walter E. Loomis, Botany , 4th ed. 
(New York: Holt, Rinehart and Winston, 1967), p. 86. 
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absorbed as light reaches the palisade layer. Green light is 
absorbed as well as red but in lesser amounts; for example, 
Wilson and Loomis report the average leaf absorption of light 
in the green-yellow region (.5 to .6 micrometers) is 68 percent 
while that in the orange-red (.6 to .7 micrometers) is 78 
percent. Thus, in the MSS data, a decrease in reflection 
will occur in Band 2 with respect to Bands 1, 3, and 4 in 
pixels dominated by actively photosynthesizing vegetation. 

As a plant approaches senescence, the chlorophyll— a deteri- 
orates and absorption decreases. ^ 

Because leaves reflect increasing amounts of near- 
infrared light as they mature and because the pigmented parts 
of the mesophyll reflect decreasing amounts of visible light 
as chlorophyll increases, the amount of reflectance in Bands 
2 and 4 would be a good indicator of vigorous vegetative 
cover. In this study, the ratio of Band 4 (second near- 
infrared) over Band 2 (red) is used as an indicator. Since 
an increase in Band 4 is associated with a decrease in Band 
2, a ratio of the two bands is more sensitive to vegetation 
changes than, for example, one band alone. Similar ratios 


Wilson and Loomis, Botany , p. 93. 

2 

For examples of reflectance spectra of plants in vivo 
at all stages of growth see S. G. Ungar, W. Collins, J7 Coiner, 
et al., preparers. Atlas of Selected Crop Spectra, Imperial 
Valley, California (New York: NASA Institute for Space Studies, 
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have been used previously in remote sensing work. 1 

Plotting the ratio of Band 4 to Band 2 on the y-axis 
against time on the x-axis for a specific 595-hectare ground 
area, the quad, yields a curve for each of the four years of 
the study, similar to the diagram below. 



It is hypothesized that a change in the shape of a plot, such 
as a shift in the peaks, from one year to the next represents 
a change in agricultural system because of the different 

temporal growth patterns of vegetable and field crops. The 

♦ 

year 1976 is particulary important in the Seabrook area 
because of the land use change described earlier; a sharp 
departure from the curves of the three previous years would 
indicate a rapid adjustment to the Seabrook Farms' decision to 
cease farming and processing. 


_ . , _ ac J ain ' Kristof and Baumgardner, "Multispectral 

Soxl Patterns, and also Eugene L. Maxwell, "Multivariate 
System Analysis of Multispectral Imagery," Photogram me trie 

and Remote Sensing 42 (September 1976): 1173-1186. 


Although measurements of reflected energy are made at 
approximately the same time Of day, the sun elevation at a 
specific latitude at a given hour changes throughout the year. 
To account for this variation, a sun elevation correction 
factor was introduced which adjusts all measurements to a 
solar angle of 90° through the equation 

sun elevation _ 1 

correction factor ~ ” ' “ ; ~ 

sxn (sun elevation) 

In addition, the use of the Band 4/Band 2 ratio reduces 
solar angle differences as well as signal fluctuations and 
differences attributable to changing atmospheric conditions 
by using a relative measure, a ratio. 

Description of Method 

The analysis employs digital LANDSAT data and relies 
upon, as support materials, LANDSAT imagery, aerial photo- 
graphy, United States Geological Survey maps, Seabrook Farms 
Corporation maps, and ground-collected information. 

Although LANDSAT data are available for passes made 
every eighteen days from July 1972 to January 1975 and every 
nine days from January 1975 through most of 1976, not all of 
it is useable „ Weather conditions such as thick haze or cloud 
cover over the study area make certain passes unuseable for 
analysis. Other technical problems occasionally preclude the 
transformation of data from certain passes into computer 
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compatible tapes. Out of 119 possible coverages, 26 acquisi- 
tions (or 22 percent) were found to be useable in this study. 
All 26 scenes are preprocessed for geometric correction and 
rectification to a Universal Transverse Mercator grid. 

Imagery is used initially to choose suitable passes and 
to geographically identify the corresponding digital data. A 
block of data over a geographic area of approximately j 4,225 
square kilometers (13,225 square miles) is called a "scene," 
and each data tape contains information from an entire scene. 
From each scene twenty-four smaller areas of 595 hectares 
(1400 acres) are used as the basic unit of analysis. It was 
felt that this unit was suitable for the long-run objective 
of the study — to simulate global monitoring. As well, the 
unit conformed to certain technical restraints in the computer 
programs which are used in support of the project. Each set 
of twenty-four areas , called "quads," is registered spatially 
from date to date. Therefore the study concerns the change 
in reflectance of twenty-four identical areas over the course 
of four years — 1973 through 1976. Map 3 shows the spatial 
distribution of the twenty-four quads in the study. 


■^For a description of geometric correction and recti 
fication, see Ralph Bernstein and Dallam G. Ferneyhough, Jr. , 
"Digital Image Processing," Photogranimetric Engineering and 
Remote Sensing 41 (December 1975) : 1465-1476. 

2 For a list of LANDSAT tapes used in this study, see 
Table 6 in the Appendix. 
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Map 3. Distribution and Land Use Designation 
■of the Twenty-four Quads Used in the Study. 










A classification routine developed at the Institute for 

Space Studies was used to obtain a nominal spectral response 

for each quad. 1 Normally, this routine is used to delimit 
clusters of pixels whioh have slmilar response curves . ^ 

this study, however, the object was to delimit one single 
cluster for each quad into which approximately 95 percent of 
the pixels fall, regardless of similarities. This percentage 
was chosen to filter out highly reflective, exotic pixels 
which might influence the quad's spectral response dispropor- 
tionately. The classification routine is applied to twenty- 

four quads from each of twenty-six passes to obtain 624 charac 
teristic spectral responses in four bands. 

For each pass over each quad, the reflectance measured 
in Bands 4 and 2 is ratioed, and this ratio is then plotted 
against time on the y-axis, with time on the X-axis. 


es ^ bli ?hing“L C nSmiifrvaiS2s r ?o“eLrclSLr e 2e d d 0f 

Algorithm, " ^in^reparation, " A £ l °S 

Earth Resources SaL InJiiS?^^ ? irector - 

New York, New York 10025 lnstltu te for Space Studies, 


i-wroT i B ji 
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RESULTS 


The results of this study are shown graphically in 


Figures 6a through 6x, which display the Band 4/Band 2 ratios 


over four years for each quad. The exact date of the LANDSAT 


pass whose data was used to determine each ratio is found either 


above or below the point corresponding to its ratio . The Band 4/ 
L . n d 2 ratio performs as was expected in that the results show 
that in all quads there is a general correspondence between the 


variation in the ratio with time and the variation in the ground 


truth and reference data. There is a depression in the ratio m 


the late Fall and late Winter-early Spring months and 


during the major planting, early growth, and harvesting 
periods due to lack of vigorous plant cover. The timing of these 
periods is documented by personal observations during the season 
of 1976, discussions with farmers, a four-year file of general 


cropping conditions in New Jersey compiled from the Weekly W e a - 


ther and Crop Bulletin, 1 the crop cycle calendars shown m Figures 


4 and 5, and information obtained from Seabrook Farms. In addi- 


tion, the two agricultural systems under study, vegetables and 
field crops, performed differently in their ratio plots. The 
twenty-four graphs follow, after which the results are described 


more fully. 


X U S., Department of Commerce, environmental Data Service, 
and Department of Agriculture, Statistical Reporting Service^ 
Weekl y. Weather and Crop Bulletin (Washington, D.C.: 1973 through 

1976 ) . 
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Figure 6f. Quad 6. 
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Figure 6x. Quad 24. 
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Comparison of Quads to Each Other 

Though the discrepancy in the dates of data acquisition 
from year to year makes statistical comparisons of the four 
curves of a quad with each other difficult, it is possible to 
compare statistically the curves of "vegetable" quads with 
those of "field crop" quads. Two such quads are shown in 
Figure 7. These quads, numbers 2 and 23, were chosen to 
represent field crops and vegetables, respectively, from per- 
sonally collected ground truth information and maps from 
Seabrook Farms. Figure 7 overlays the two types of systems so 
that their differences are readily observable. Table 7 in the 
Appendix lists the exact differences between the ratios; 
generally, they are most distinct in July and August, with an 
occasional distinction in April, May, and June. 

The calculation of the Band 4/Band 2 ratio was extended 
to eighteen other quads labelled either vegetable or field 
crop based on the information at hand. Using the nonparametric 
Mann-Whitney "U" test, the two groups were tested for differences. 1 
The means of the two groups were significantly different on 
data from nine of the twenty-six acquisitions, with at least 
a 95 percent level of confidence. Seven of these nine passes 
were during July and August. Again the summer months show a 
significant divergence between agricultural systems. The "U" 
scores for each date are listed in the Appendix in Table 8 . 

1 Sidney Siegel, Nonparametric Statistics (New York: 
McGraw-Hill Book Company, 1956 ), pp. 116 - 127 . 
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The exact differences betvzeen the vegetable and 
corn/soy ratios are listed in the Appendix in Table A2. 
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Although the ratios tend to show differences at a specific 
time of year, the reason proposed for these differences has 
more to do with land utilization practices in space rather 
than in time. The curves diverge within the summer months 
when both field crops and vegetables are at a maximum, perhaps 
the very time when convergence would be expected if one con- 
sidered crop cycles alone. Field crops tend to be in full 
canopy during July and August whereas vegetables tend to have 
less canopy cover with rows being separated by strips of soil. 
The large amount of exposed soil decreases the Band 4/Band 2 
ratio, and this is evident in the figure. 

Curve separations or similarities may also be related to 
meteorological events which immediately precede data acquisition. 
For example, the LANDSAT pass during the Summer of 1975 
(August 2) occurred ten to fifteen days after a series of 
severe storms destroyed crops in the southern New Jersey region. 
On Figure 7, this event is reflected in the similarity of the 
two plots of 1975 during the Summer, when they would be expected 
to diverge. 

The ratios for the year 1976 do not behave differently 
from those of the preceding years. Two explanations are pos- 
sible, namely (1) the transition from vegetables to corn and 
soybeans, complete in 1977, may not have been far enough ad- 
vanced in 1976 to be detected; or (2) residual factors of 
vegetable farming practices may still be impacting former 
Seabrook Farms land in 1976 to some extent. Supporting this 
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second possibility, the former Seabrook Farms' Farm Manager, 

Mr. Val Perri, is of the opinion that the strips formerly used 
in rotation of vegetable farming will be evident in the field 
cropping system for a few years to come.^ 

To test the separability of a co-existing natural vege- 
tation component in the landscape, an analysis of variance was 
performed using the nonparametric Kruskal-Wallis test on three 

types of vegetative systems — vegetable, field crop, and riverine/ 
2 

forest. Quads containing riverine and forest vegetation also 
contained agricultural land, but they appeared to group together 
in the data. Observe the shapes of riverine/f orest quads 5, 

9, 13, 14, 17, and 22; these quads have the widest range of 
ratio values and form simple plots with few oscillations. The 
Kruskal-Wallis statistic revealed that the three sample groups 
were most probably (99 percent confidence for June 1973 and 
95 percent for June 1976) taken from different populations. 


Two-Year Cycle 


A two-year cycle is evident in some quads. In twelve 
quads, the plots for 1974 and 1976 were similar. 2 * Table 7 



^Val Perri, Seabrook Foods, Incorporated, Warehouse 
Manager, personal communication, August 4, 1977. 

2 Siegel, pp. 184-193. 

3 ... 

The similarity of the plots was determined by the agree- 
ment of two persons who visually compared each of the years 
within a quad to each other. When both persons felt that 
the plots of two years were alike, they were designated "simi- 
lar" and entered as such in Table 7 . 



TABLE 5 


Quad 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


PAIRS OF YEARS IN EACH QUAD 
WHICH HAVE SIMILAR PLOTS 1 


73-74 


73-75 


73-76 


74-75 


X 


X 

X 


X 

X 


X 

X 


X 


74-76 


X 

X 

X 

X 


X 

X 


X 

X 


X 


W "similarity" of the plots was determined by the 
agreement S two'peLons based on visua! inspects. 
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lists all quads and their year-to-year similarities. This 
pattern may be explained by the wheat-soybean cropping system 
which many farmers practice. In this system, wheat is planted 
in the Autumn and harvested in late June or July. A late crop 
of soybeans is then planted and harvested in late October. 

By this time the growing season is too short to drill in a new 
crop of winter wheat and the land remains clean-tilled. The 
following Spring an early crop of soybeans is planted, and this 
is harvested in time for a fall planting of winter wheat, two 
years after the first winter wheat crop. Not all wheat-soybean 
fields in a quad will be in the same stage of the cycle at any 
one time, but the effect is still present unless exactly half 
of those fields are in one stage and half in the other. Theore 
tically, a half-and-half situation would cancel the effects of 
both stages of the cycle and the cycle would not be evident. 
However, it is unlikely that this is the case, and thus the 
■j- W o— year cycle is exhibited. 
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DISCUSSION 

The method used in this experiment to monitor agricultural 
change and its approximate location is based on a sound prin- 
ciple; that is, the identification of agricultural systems by 
tracing the periods of peak living biomass in a fixed area over 
the course of a year. Some problems exist with the method as 
it was used in this experiment which prevented optimal results 
and which should be noted when extending the method to other 
areas. Some of these problems relate to lack of ducumentation 
on certain cropping practices when the method was devised; and 
others are instrument-dependent, relating to the temporal reso- 
lution of the satellite, the impenetrability of cloud cover, 
and the occasional failures of telemetry and data processing. 


During the winter months the Seabrook land was planted 
to a cover crop of rye or wheat, which was cut before the crop 
reached maturity. Then the land was turned and vegetables 
were planted. Data collected in short time ihtervals — nine days 
should be suf f icient--would exhibit the abrupt shift from 
fairly thick green vegetation to bare soil. But if the data 
collection platform missed the cutting and cultivating period 
(because of atmospheric or mechanical problems) , the results 
could be misleading to the analyst. 

If winter months are used in the separation of agricul- 
tural systems, cover crops may be easily confused with winter 
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grains. In addition, there is some fall planting of vegetables 
which might be confused in the October reading with cover crops 
and winter grains if the grains were planted early. 

The growth patterns of corn and soybeans are similar and 
are here considered to be part of the same agricultural system. 
Their growth patterns are distinct from other crops: planting 

in May and June; three to five leaves showing through late June; 
luxuriant green growth in July and early August; yellowing 
through September and October. A late crop of soybeans planted 
after wheat might confuse the signal in terms of timing although 
not in terms of sequence. In September and October, however, 
it should be fairly easy to watch the ripening of corn and soy- 
beans with this technique. At that time, either fall vegetables 
(for example, fall mustard and turnip greens, squash, and fall 
spinach) or cover crops will be well along and would not be 
confused with the corn and soybeans. 

Orchards are a source of confusion since between each 
row of trees are about twenty feet of nonorc*.ard growth, usually 
soybeans or hay. If the interarbor spaces are planted with 
soybeans, they will contribute to the assumption that a quad 
is predominantly devoted to field crops, while the fact that 
orchards make up a large part of the quad will go undetected. 

That a two-year cycle of wheat and soybeans is practiced 
in many fields means that the analyst must construct his 
models of cropping around a biennial rather than an annual 
sequence of changes in the ratio. It is not known how widely 
this practice is used and so its effect cannot be estimated. 


d-P-' 




Even discounting the biennial pattern effect, it is possible 
that the wheat/soybean curve would be significantly different 
from the vegetable curve under optimal conditions of data 
acquisition. 

In reality, the temporal resolution of the LANDSAT data 
available at the user end of the system turns out to be too 
coarse for this type of analysis. Because not all LANDSAT 
passes are useable due to weather or mechanical difficulties, 
there are wide time gaps between measurements in some years 
and a lack of concurrence between years. For instance, 
except for early/mid June, there is no time period which is 
represented in all four years. 

If one were to extend this method to other areas of the 
world where different types of agricultural systems existed, 
certain adjustments should be made. Smaller or larger field 
sizes than those in southern New Jersey might justify an 
increase or decrease in the size of the areal unit used. 

A higher degree of crop mixing might mean a broader definition 
of "agricultural system" would be in order. Cyclical charactens 
tics of the crops grown merit attention. Finally, when two 
agricultural systems are being compared, the seasonal aspects 
of the systems must be considered to determine if nonoverlapping 

periods exist. 
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CONCLUSIONS 

A change in agricultural practices in an area of southern 
New Jersey has offered the opportunity to experiment with 
agricultural systems identification methods using LANDSAT. 

The method described in this paper is a promising tool 
because it differentiates the vegetable and field crop systems 
and because the shape of the curves over time generally 
corresponds to growth patterns. Summarized, the experiment 
revealed the following conclusions: 

1) The Band 4/Band 2 ratio, when calculated for 

a 595-hectare area, differentiates vegetable from field 
cropping systems most effectively in the months of 
July and August. 

2) Areal units of the size used in the experiment 
are also sufficient to identify regions where riverine/ 
forest vegetation predominates. 

3) A biennial cycle, where it exists, must be 
incorporated into the analyst's hypothesis. 

4) In the agricultural systems investigated, 
patterns of land utilization in space rather than in 
time have the most impact on the behavior of the ratio. 
With better processing technology or microwave sensing, 

it may be possible to eliminate large gaps in the series of 
data acquisitions. Microwave sensing eliminates the problems 
imposed by cloud cover; however, different indicators sensitive 
to vegetation must of course be developed in the microwave 


portion of the electromagnetic spectrum. 

A LANDSAT pass every nine days would be sufficient for this 
type of experiment, were each pass useable. A more frequent 
data collection system would increase the probability of passing 
over a given area on a cloud-free day. Improved similarity 
between the dates of acquisition in each year would increase 
the ability to detect a change in the shape of a curve from 
one year to the next, and it may also permit the use of more 

formal analysis to quantify this difference. 

The results from the study are encouraging enough to sug- 
gest that the method employed would be sensitive to a change 
from one agricultural system to another. The methodology may be 
useful in other areas of the world undergoing agricultural 
transition with field sizes and a degree of crop mixing similar 
to that of southern New Jersey. Alternatively, the method 
is adaptable to the agricultural systems which are being 

monitored. 





TABLE 6. LANDS AT DATA USED IN THIS STUDY 


1973 


1974 


1975 


1976 


DATE 

SCENE 

IDENTIFICATION 

NUMBER 

February 13 
June 1 
July 7 
August 3.2 
August 30 
October 23 

1205-15141 

1313-15141 

1349-15134 

1385-15131 

1403-15125 

1457-15113 

April 3 
April 21 
June, 14 
July 20 
October IS 

1619-15083 

1637-15080 

1691-15063 

1727-15052 

1817-15023 

May 22 
June 9 
August 2 
October 13 
November 18 

5033-14530 

5051-14520 

5105-14485 

5177-14443 

5213-14422 

March 5 
March 23 
April 10 
April 19 
April 28 
June 12 
July 18 
August 23 
September 28 

5321-14352 

5339-14341 

5357-14330 

2453-14571 

5375-14314 

2507-14553 

2543-14545 

2579-14535 

2615-14525 
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TABLE 7 


THE RATIOS OF QUAD 2 (FIELD CROPS) 
MINUS THE RATIOS OF QUAD 23 (VEGETABLES) 


DATE 


DIFFERENCE 


1973 February 13 
June 1 
July 7 
August 12 
August 30 
October 23 


- 0.10 

-0.29 

0.88 

1.70 

1.23 

-0.24 


1974 April 3 
April 21 
June 14 
July 20 
October 18 


0.33 

- 0.01 

-0.58 

1.94 

omitted 



1975 May 22 
June 9 
August 2 
October 13 
October 31 
November 18 


0.99 

1.27 

0.11 

0.70 

0.53 

0.36 


1976 March 5 
March 23 
April 10 
April 19 
April 28 
June 12 
July 18 
August 23 
September 28 


-0.03 

0.08 

-0.18 

-0.54 

- 1.00 

-0.27 

2.32 

1.13 

-0.26 


TABLE 8 


TESTING THE DIFFERENCE BETWEEN THE MEANS OF THE 
VEGETABLE AND MELD CROP QUADS USING THE MANN-NHITNEY "U" TEST 


DATE 


"U" SIGNIFICANCE 3 

STATISTIC 


1973 February 13 
June 
July 7 
August 12 
August 30 
October 23 


27 
30 
15 

3 

3 

28 


.05 

.002 

.002 


1974 April 3 
April 21 
June 14 
July 20 
October 18 


25.5 
24 

37.5 
16 

mitted 


.10 


1975 May 22 
June 9 
August 2 
October 13 
October 31 
November 18 


18 

27.5 
8 

14 

21 

26.5 


1976 March 5 
March 2 3 
Paril 10 
April 19 
April 28 
June 12 
July 18 
August 23 
September 28 


' 22 

27.5 

34.5 
34 
21 

29.5 
5.5 
8 

28.5 


.02 

.05 


.002 

.02 


a For a two-tailed test 
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MONITORING AGRICULTURAL SYSTEMS* 


Patricia J. Ncary and Jerry C. Coiner 

Columbia University, Department of Geography, and 
t NASA, Institute for Space Studies, New York, NY 10025 

To a large extent, remote sensing of agricultural phenomena has 
concentrated on analysis of agriculture at the field level. Concern has 
been to identify crop status, crop condition, and crop distribution, all 
of which arc spatially analyzed on a field-by-field basis. This is not 
the only geographical scale at which agricultural activity can be analyzed. 

A more general level of abstraction is the agricultural system, or the 
complex of crops and other land cover that differentiate various agri- 
cultural economics. This paper reports on a methodology to assist in 
analysis of the landscape elements of agricultural systems with LANDSAT 
digital data. These methods may be particularly useful in the context of 
the developing nations to monitor agricultural change associated with 
economic" development. 

The site used to develop and test this method is part of Cumberland 
and Salem counties. New Jersey, in the northeastern United States, where 
change from a traditional commercial vegetable system to a corn-soybean- 
wheat field crop system is taking place. During the period of transition, 
both systems exist in close proximity, as farmers make individual decisions 
about what crops to grow. The method to differentiate these two systems 
relies upon an indicator which measures vegetative cover on a given ground unit 
during the course of four years. The unit employed is approximately 595 

*This research was supported by NASA Grant, NSG-5080. 



hectares and way encompass up to 40 fields. Units representing vegetable 
cropping were chosen from corporate farm land specializing in vegetables, 
while units representing field crops were chosen from private landholdings 
planted to corn, soybeans, and wheat nearby. 

The method of analysis tested in this experiment is based on the 
theory that, since agricultural systems differ in their crop calendars 
and cropping practices, the amount of vegetative cover in a given area 
over the course of a year will depend on the cropping system employed. 

The method involves tracing the periods of photosynthetic activity in a 
fixed area during a four-year period to determine if land devoted to 
field crops can be, differentiated from that planted to vegetables at the 
595-hectare level of aggregation employed here. 

The measure used to determine the presence of green vegetation on an 

unit of land is a ratio of LANDSAT Band 4 to Band 2 (Band 4 ranging from .8 

to 1.1 micrometers, Band 2 from .6 to .7 micrometers) in terms of milliwatts 

2 » 

per steradian-centimeter . This ratio is particularly sensitive to pre- 
senescent vegetation since Band 4 reflectance increases with the development 
of the spongy mesophyll section of leaves while Band 2 reflectance is de- 
creased due to chlorophyll absorption. The plotting of this ratio over time 



Figure 1. Example of the change in the Band 4 /Band 2 ratio over time. 
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The results from two units, one predominantly vegetables, the other 

field crops, are shown in Figure 2. The shapes of the curves differ from 

year to year because of a lack of concurrence in satellite acquisition times 

' for usable data. The areal units can be compared to each other within 

each year, however. The curves diverge within the summer months when 

* 

both field crops and vegetables are growing. Field crops tend to be in 
full canopy during July and August whereas vegetables tend to have less 
canopy cover with rows being separated by strips of soil. The large amount 
of exposed soil decreases the Band 4/Band 2 ratio, and this is evident 
in the figure. Curve separation may also be related to meteorological 
events which have immediately preceded data acquisition. For example, the 
LANDSAT pass during the summer of 1975 occurred ten to fifteen days after 
a series of severe storms destroyed crops in the southern New Jersey region. 
On Figure 2, this event is reflected in the similarity of the two plots of 
1975 during the summer, when they would be expected to diverge. 

The calculation of the Band A /Band 2 ratio was extended to twenty-two 
other 595-hectare units. Each unit was assigned to a vegetable or field 
crop category. Using the nonparametric Mann-Whitney "U" test, the means 
of the two groups were significantly different on data from nine of the 

J 

twenty-five acquisitions, with a 98 percent level of confidence. Seven of 
these nine passes were during July and August. Again the summer months show 
a significant divergence between agricultural systems. 

To test the separability of a co-occurring natural vegetation component 
in the landscape, an analysis of variance was performed using the nonpara- 
metric Kr u ska 1-Wal 1 is test on three types of vegetative systems— vegetable, 

v 

field crop, and rlverine/f orest. In both June 1973 and June 1976, the 






statistic revealed that the three sample groups were most probably (99 per- 
cent confidence for June 1973 and 95 percent for June 1976) taken from 

different populations. 

As a preliminary test of an aggregated-unit approach to monitoring 
agricultural, system, the results are encouraging. The LANDSAT-derived 
radiometric indicator, even at this level of aggregation, appears 
to afford a significant means to differentiate agricultural systems from 
each other as well as from natural vegetation elements. Future experiments 
will be performed using sites in both temperate and tropical regions, 
to extend the utility of the indicator in monitoring agricultural systems. 
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Given a need to collect serial data for regional 
planning, the purpose of this paper is to describe a 
method which uses both remote sensing and information 
systems for collection, spatial analysis, and display 
of periodic data. Findings presented hero are from 
an ai rex *.f t/Landsat based change detection study con- 
ducted for the period 194C-1972 on Marinduque Province, 
Republic of the Philippines. The experimental informa- 
tion system, for the Island allows collection and 
maintenance of spatial relationships inherent in the 
landscape through generalization of a preselected 
! number of data elements into a uniform grid cell 

(twenty-f j.ve hectares). Storage is achieved by array- 
, ing the data elements as records, based on the grid 

1 cell's location within the predefined Universal Trans- 

• verse Mercator (UTM) geographic reference system. 

. Either aggregated data, categorized and combined into 

i tables, or spatially displayed data in the form of 

! • computer maps, can be recalled for analysis and plan- 

' nibg support. Programming allows presentation of 

; dis crete symbology or choropleth categories in computer 

i generated maps . 

: Each of the 4,008 twenty-five hectare cells repre- 

| senting Marinduque were studied for changes in and 

between various development variables. These changes 
j were measured and tested using nonparametric statis- 

' tics to determine the effect of specific land cover 

! . changes on other infrastructure components and proxi- 

i mate land use. Since the grid was spatially registered 

through time, change occurring within the cell could be 
j ’ positionally related, allowing assessment of the 

| change's impact. 

' Techniques have been developed for more continuous 

] updating of the Marinduque data base with Landsat data. 

| The primary effort has been to extract for each 

1 twenty-five hectare cell an indicator of change as 

! derived from UTM geometrically corrected and aggregated 

j handset data. Anomalous or unique changes over time 

! can thus be detected and corroborated either with 

, ground , reconnaissance or high resolution aircraft 

imagery. 


The experimental system used here demonstrates that 
it is possible to retain in one information system 
historical and current physical aspects of developing 
regions . In support of the planner's requirement for 
periodic observation, the study reveals both trends 
and seasonal and other periodic cycles in regional, 
settlement, and environmental patterns which, x£ 
unaccounted for by planners, would impact negatively 
on development projects. If properly used, these data 
Drovide better regional models than static base line 
approaches, thus markedly improving present development 
planning capabilities. 


1. INTRODUCTION 


are short, of unknown rcl-u -l ty is compounded by a lack of timely 

United .lotions, Lerh-inisms which make much of the data useless for 

collection nnct analysis mechanisms ‘ t-houah data deficiencies are 

Mho.' than is 

only One aspect o_ the i mhci'miHallv to the failures of development 

aouhh Ihijt the, can coht.'.oot. .ob .973) . This lack »« 

...proachos, suoh OS thoso aro' anavaiiablo 

thing that is uncommon even in developed countries^ O^P sensi sources to 
the problem data in development projects _ WAS , 

; s sensing, this paper presents «n experimental 

geographic information System which supports the various aspects of the plan 
ning process. 

Before considering the experimental system one of J £ a “mSto 

tions about data from remote sensing, now y . l’eneed (Landini and 

sensing geographic information * ^ ed ? 9 ^ ) ^ C ?he assumption is that 

remotc^scn Ling” data acquired' at a single time and used to 
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assumption is subtly misleading , and it rc •' tha static data base 
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. *? St th ® utility of remote sensing and a geographic information system 

in a development context, a study site was. chosen in the Philippines. The site 
Marinduque Island, Mannduque Province (see Figures 1 and 2) is a tvoical 
underdeveloped area prevalent in the tropical world. It displays the statis- 
tical symptoms of underdevelopment from the structure of its population— 
predominantly rural the uses of its resources— generally extractive and 
unrenewable— and the basis of its economic activity— largely subsistence agri- 
culture and lacking alternatives (Coiner, 1975) . 9 

2. INFORMATION SYSTEM 

The number of features that can be interpreted from aerial photoqraohs or 
extracted from Landsat data or digitized from more conventional data sources 
(maps and census series) require design of an information system to manage the 
large anc. complex cata sources. Information systems delimit an organized 
approach to the specification, collection, storage, analysis, and display of 
data. When spatial referencing is contained within the data file structure 
gt^ir') , r te i'Qif\ kn0 c n ? S a ‘J e °9raphic information system (Tomlinson, 1972; 

/h 1974 ii f UCh a ®y st ? ln has been designed to support this research. It 
allows the collection and maintenance of spatial relationships inherent in the 
landscape through generalization of a preselected number of data elements into 
a uniform size cell 25 ha). Storage is achieved by arraying the data elements 

Serenfnv ° n ° r - d Cell ’ S location "ithin a Universal Transverse 

^ . 1 l 9 coreferencing system. Either aggregated data, based on cate- 

9 ?* on and combination of the data fields, or spatially displayed data, in 
the form of computer maps# can be recalled for analysis (Figure 2) . 

Computer cartography for this system is based on the use of an IBM 370/145 
and a CALCOMP drum-type plotter. Programming allows p-esontat.TonCf discrete 
synOiology or choropleth categories, and display of either symbol o? patterns 
can be accomplished in multiple colors. ^ 

of A-ln® J???® !? aSe W<1 ! constructed using a 500m x 500m grid as the basic unit 
of data collection and storage. For each of the grid cells, data were inter- 
format in^able 1 !* 1 SCale aerial P h °t°grapby on a cell by cell basis using the 

fh , ( . ” a ^ ndUqU - ^° n ’? U . t f r files . arc different from traditional base line data in 
consist of three aerial photographic coverages (1948, 1967, 1972) 
narts a , tl f SCL ‘ ics ; b ?ch of the three data sets were divided into two 

land 1 co l VC:r classification, and 2) specific data related to develop- 
ment, particularly human settlement infrastructure. Data elements were H 

road? nC or t ^^ l T catc * ories ^ s P^tial type, point (e.g., house), line (e.g., 
roacij # or area (e.g.# mangrove swamp), 

Tc1 -.^ and C ° Ve ? ? at fJ° r i es interpreted from aerial photographs of Marinduque 
Island were originally developed as part of an effort* to support agrarian 9 

Bruce™ igrr Ug v° Ut ^ he Ph i 1 i L l Ppines ' Particularly Central Luzon (Veracion and 
cSssilM^t Veraclon «nd Bruce, 1967). Because of its original purpose, the 
gorics modification and supplementary land cover cate- 

b f Ma f indu g ue s rte. Twenty-one land cover categories were 
numeric oLp S *-o rCa /i e ^ tenSiV ?- 0n Marinduc ^ ue (Table II) and were assigned a 
Wic L ! h al ?, ln creating computer compatible files. The classification 

nenernl A three-digit computer code was as follows: hundredths digit — 

'Si c ° v ® r u Vy pe; ^nths and lower digits— specif ic land cover. This 

Had simil-,r i Tn t< i d hlelarchlc al classification structure. To mark areas that 
had similar land covers, a photo symbol code was applied to the photographs, a 

pretation *7, ? ffV ^° ff WaS 3 separate operation from P Uoto inter-' 
recreant' Lac l 1 f f !" he 4 '°° 8 twenty-five hectare cells, which in combination 
virtual map of Marinduque Island in the data base, then could be 
ln a u number of w «y s - Since the grid was spatially 
Slated and evaluated? 6 ' 9 * ° CCUrrin 9 within cell could be positionally 
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UH 1 file structure of the data elements interpreted from the aot-ini 
photography provided the framework for incorporating Landsat data into the 
information system. Conceptually, concatination of Landsat data into in mri-ii 
photographic data base is accomplished by identifvina an 

the Landsat data which represents the same ground^rea (here?" 25 £ ha cell ) W in Hln 
both sources. However, the Landsat data first must be reprojocLd into the V, 
map referencing system. For this study, tho renroieet-inn 'i UT 

geometrically correcting the digital data durin? preprocessing (Eernstoi^aSd 
Fernoyhough, 1975, pp HG5-76). Using the Institute Studies f acu- 

ities, a maximum of S3 geometrically corrected pixels (7 x 9 pixel mar riv? 1 
comprised a 25 ha cell, if partial pixels were included, and iftheviere 
excluded, a minimum of -IS pixels (G x 8 pixel matrix). Figure 3 is a nrmlhi,. 
representation of the congrucncing process. The air phote on the left in the 
figure was overlaycd with a 500m UTM grid taken from a 1:50,000 base map On 
the right, the Landsat HS5 band 7 (geometrically corrected) CRT imaae also w.c 
overlaycd with the same 500m reference grid. The land cover attributes o? each 
coll in the photograph wero extracted by qualitative interpretation while for 
the Landsat data, the same area was described numerically in terms of radiance. 

-*)«/£»£ SlSrt" ) l,e 

ss“. itc ”? in f° •'> chi».T?;aic2« 

(/G.20m X GO . 9Gm) , change duo to locational error was minimised Also the 
selection of a 5 x 7 pixel matrix within each 25 ha unit minimised certaS 
typos of boundary effects. Although choice of a Landsat radiance descriptor as 
a change indicator for each coll remains a matter for conHnn?n« 1 f. 

data in each MSS band -ore described by three statistical n i- ! m ' th ° 

standard deviation , and percent standard deviation, other properties of ^t-he' 0 ^ ' 
matrix may be equally or more appropriate. properties oi the 

Since the airphoto file and the Landsat file were linked bv loraHrmai 
properties, they could bo used to support each other. For instance the 1972 
aerial photographic data base could bo used to identify colls within the 19?2 
Lanco.it file which had a homogeneous land cover (c.g,, paddy rice, coconut 

^Ssod*te 0 |^tibliSh r ?,«rt?;{ Prty “^“r* 0rCB e r *9* owth) - TheU ceils then could 
,7°, ®,V" L ° Gotnblioh opcctral properties of various land covers. Cells which 

. d .diffeient. secondary and tertiary land covers could bo included uo studv the 
point data^el emont<'°earri‘ert **1®* "’^ th l uixt r d land covers. Further, line and 

=?*£ sssfeo tdTitiiz * cor 

their aggregate statistical impact on a 25 ha area. J univs) trough 

In short, these types of analyses permit use of two da«a sources in 

increasing the utility of both. Aerial photography provides historic 
cycles "and *"t inie 1 v''da t a on'lmtri Provides data to study short-term environmental 
activities. ly d 1 d ‘ changcs whJch »ay relate to developmental 

| ANALYSIS OF LAND COVER CHANGE FROM AERIAL PHOTOGRAPHS (19-18-1972) 

tion change can be cyclic, e.g., tho encroachment of parang vegeta- 

tion on coconut orchards, the subsequent removal of the pa van a and it's 

innovation** "nth pla ’/ l ' ilKJU ' or noncyclie, e.g., introduction of an 

lLnd cover * n i- * opme V^ P«>Dect, which manifests itself as a new 

tion ^ c ?*?. of cyclic change, the periodic expansion and contrac- 

rnen ai or lnman b ? tw ?, on lfmd c ° vorn * generated by either environ- 

activities. In tho case, of developmental change, the impact of 

expense "of the' old. P ° thc Oxpansion of * land cover at the 

Since both cyclic and developmental (noncyclie) land cover change may take 

firm G o? t-Se in a «^cific place requires some 

inis „ 7 C1 x p- Analyses to measure tho impact of the two land cover 

cover change ^ 'the 'techniones 08 ! an y\ y f of! ha y? not boon used widely to study land 
ci change, the techniques should be applicable. The innovative component of 
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land cover change would be measured by the linear trend in the data, while the 
cyclic would be represented by the periodicity of the regression error. 

A third type of change within the land cover which may reflect development 
activity is land cover fragmentation. As land covers are added to a small 
fixed areal unit, such as the 25 ha cell in those files, it is recorded as an 
increase in the number of land cover types occurring in a given cell, and 
fragmentation results from such increases. For example, in a three category 
cys tent, each cell could contain up to three cover types with a maximum of 
six (3!) possible land cover combinations. When just one new land cover is 
introduced by & development prefect, the number of possible combinations 
increases to twenty-four (4*). This leads to the hypothesis, testable with 
files used in this study, that development projects will cause land cover 
fragmentation by the introduction of new cover types into the landscape. 

between 1948 and 1972, the number of land cover types on Marinduque 
increased from 15 to 23, an increase of 8 new land covers during the twenty- 
four year period. The greatest increase in land cover categories was in the 
areally secondary class, where the number changed from 13 typos in 1948 to 22 
in 1.9 12 . 1 ho areally dominant class changed least, with only 6 now cover types 

added (Table III) . Land cover types which were identified as new in 1967 and 
19/2 are indicative of the physical development that had occurred. These 
include now transport modes (airports), new agricultural activities (fish ponds 
and sugar cane) and now industrial and mining activities (strip mininq, indus- 
trial sites). 

A land cover map (Figure 4) of Marinduque provides an assessment of the 
spatial distributions of dominant land covers by 25 ha cell for 1972. In 
addition to static land cover distributions, land cover change can be measured 
in terms of the number of cells dominated by a given land cover nr.d the varia- 
tions in this dominance from one time to another (Table IV). A significant 
difference was found in the dominant land cover distribution for the 1948 to 
1972 period. This difference was also evident in the north central-northeastern 
portion of the Island for the 1948-1972 period, although the Chi-square value 
was considerably smaller. Testing of the 1943 to 1967 dominant land cover 
distribution also revealed a significant difference. In the shorter 1967-1972 
period, tlie distribution of cells dominated by a given land cover type was not 
significantly different. 

The Chi-square test results revealed the degree to which each of the 
dominant land covers had undergone change. Clearly, the tests showed that the 
land cover was gradually changing, with predominance shifting from grass and 
brushland (1948) to coconut (1972). 

The above analyses concentrated on the dominant land cover in each cell 
and disregarded secondary and tertiary land covers. However, evidence in 
Table III indicates that more now land cover types were introduced as secondary 
or tertiary land covers than as dominant ones, tncrcforc, more change could be 
expected in these categories. 

Changes in land cover can be determined by constructing an array of 
possible combinations (dominant land cover, secondary, tertiary or new) and by 
inspecting each cell for each time period and assigning it to a particular 
change state, i.e., no chango, shift in predominance of existing land covers, 
addition of new land covers. Analysis shows that the majority (85 percent, 
1948-1972) of cells underwent some type of land cover change. The most common 
change , which occurred in 17.3 percent of the cells', was the introduction of a 
new tertiary land cover. Complete change of all cover types in a specific coll 
was rare, occurring in only .5 percent of the colls in the 1948-1972 period. 
Cases where existing land covers in a cell shifted in relative importance 
totaled 35 percent of the cells. Spatial patterns of change shown in the map 
(Figure 5) reveal that areas where new land covers (referred to as "land uses” 
on the map) have been introduced ring the Island's seaward-facing lowlands, 
center on the Island's all-weather road, and include the upland area where the 
copper mine is located. Cells that have undergone land cover shifts from one 
degree of dominance to another tend to be located in the upland interior of the 
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Island whc-re the human population is less dense and where swiddoning is prac- 
ticed. No-change colls are confined to two situations: (1) core, areas of 

human dominance/ around the poblaaiono (county scats) on the heavily settled 
lowlands, and (2) inaccessible areas, dominated by natural vegetation, which 
hinder human entry with physical barriers, c.g., steep slopes, swamps, etc. 

When the three most prevalent cover types were analyzed together, a 
upecific series of land cover combinations were formed. To study the expansion 
and contraction of these triformed land cover units, representing cyclic 
change, those occupying more than 28 colls were identified for the north 
central— northeastern area of the Island* In 1940, there wore 11 triformed land 
covers. Of these 11, 6 showed a cyclic trend in the twenty-four year study 
period. Four additional triformed land covers expanded to occupy the minimum 
28 cells by 1967. Of these, 2 displayed cyclic trends. One additional tri- 
formed unit exceeded the minimum area in 1972 and was of a cyclic nature. 

Thus, expansion and contraction in a cyclic pattern were clearly evident in 9 
of the 16 triformed land covers. 

Although the limited number of observation points (three) in this study 
minimizes the utility of time series analysis, a simple linear regression 
model and a calculation of a Pearson product-moment correlation (r*) were 
applied to the data. Table V shows the results of linear regression on the 16 
triformed land cover units previously discussed. 

Results of the regression and correlation analyses indicate that tri- 
formed land cover units which have positive slopes (increasing area) and high 
r 2 usually contain at least one category of paddy and cogon, showing the expan- 
sion of tlie messicol land covers at the expense of semi-natural covers. Tri- 
formed units containing dry crop and oogon in combination show a negative 
slope (decreasing area) and/or relatively low r 2 values, expressing greater 
cyclic change and a total decrease in slash and burn (dryland crop) on the 
Island. Also, natural and semi-natural triformed land cover units tended to 
show negative slopes and relatively high r 2 values, emphasizing the constant 
decrease in area not directly affected by man. 

It has been hypothesized that as development causes the introduction of 
nev; land cover categories, these new categories will cause greater fragmenta- 
tion in the landscape, creating a larger number of new land cover combinations. 
Inspection of the Marinduque land cover files reveals this to be the case. The 
number of land cover combinations increased from 2B4 in 1948 to 373 m 1972, or 
31 percent. For the entire Island, combinations including urban and mining, 
fish ponds, riverine, coconut, and forest exceeded the average rate of increase, 
while the number of land cover combinations with dryland crops actually 
decreased. 

Increases in combinations referred to earlier may be the outcome of two 
‘land use systems (traditional and modern) coexisting on the Island. As the 
data illustrate, cells which have modernized land covers, particularly urban 
and mining, usually do not include traditional covers. Evidence to support the 
parallel, but separate, coexistence of the two systems includes the observation 
that the majority of new land covers held a tertiary position, meaning that new 
land covers occupied the least amount of area. A second point is that by plot 
ting the triformed land covers as connectivity graphs (Figure 6), the land 
covers introduced since 1948 tended to be isolated whether or not they were 
subsequently of tertiary, secondary or primary spatial importance. This is 
illustrated in Figure 6 by the new land cover "fish ponds" (type 109), which 
was poorly connected to other land covers in the graphs, even though over time, 
fish ponds increased spatially from a tertiary to a secondary position in some 
of the cells. 

4 . CHANGE DETECTION WITH LANDSAT 

Analysis of development related land cover change through use of aerial 
photography is dependent on irregular historic coverages, often resulting from 
high priority development demands which require an expensive aerial photo- 
graphic mission. These constraints on the acquisition of aerial photography 
limit its utility as a source of data for change detection studies. For this 
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reason, recent research has been oriented to incorporating handset data into a 
land cover change detection scheme. 

Usable Lands at data were first acquired for Marinduque Island on 12 Oct 
19172 (handset scone E-1081-01462) . The handset acquisition was nearly simul- 
taneous with the most recent airphoto coverage used in these studies and, thus, 
the logical basis for initial methodological studies to establish change 
indicators (see Section 2.). The first steps involved assessing the indica- 
tor's ability to discriminate tropical land cover types within the same scene 
(Tnb)c II). Using classes established for the photo interpretations, 25 ha 
cells which contained a single (dominant only, Table I) land cover were iso- 
lated in the base file. Twenty-five ha matrixes then were identified in the 
digital handset data, which corresponded positionally to the single land 
cover's location in the base file. Means, standard deviations, and percent 
standard deviations were extracted to use as pure land cover type change 
indicators. 

Table VI presents the pure typo indicators for the 1.2 Oct 1972 coverage. 
Inspection of the table reveals that the three indicators separate the land 
cqvers successfully, studies now in progress will expand the number of land 
covers considered, and extend the analysis to two more handset scenes, 
E-1153-01465 and E-2418-01253. 

Scene K-1153-01465, acquired 23 Dec 1972 is px-csently being processed 
to provide a similar data set to that available Cor the 12 Oct 1972 scene. 

This will permit introduction of short-term change data into the analysis, 
making possible study of. change indicator performance on a seasonal basis. 
Additionally, the intention is to incorporate long-term data from handsat 2 
scene E-2418-01253 acquired 15 Mar 197G. With those throe scenes, it will be 
possible to gain insight into both short and long-term variation in land cover 
reflectance. 

handset data have been limited in this study by cloud cover over the area, 
reducing the number of pure cover colls available for analysis, and by the 
nature be foreign data acquisition, which is highly concentrated rather than 
spread throughout the satellite's operational life. These limitations reduce 
the seasonal similarity of the data and introduce variations which make 
year-to-year comparison difficult. 

Further research, with mixed airphoto/Lnndoat files created for this 
project, will bo directed at definition of a change detection system. Present, 
work has shown that the aggregated handsat data, structured into change indi- 
cators, can be used to distinguish major tropical land cover types. Research 
indicators such as these, or others under development, should allow satellite 
systems to play an important role in monitoring development ns expressed in 
land cover change. 
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Figure 2. MARINDUQUE ISLAND, MARINDUQUE PROVINCE 
(Computer base map showing municipal boundaries, place and locational names) 
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Figure 3. UTM GRID OVERLAYED ON AIRPHOTO AND LANDSAT IMAGE 

(Distortion of grid due to displacement in airphoto and geometry of 
CRT bandcat band 7 image) 
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Figure 5. PATTERN OK LAND COVER CHANGE 1948 TO 1972 
(Differences in land cover plotted as a three category pattern map) 



Dominant 

Secondary 

Tertiary 

LAND COVER 


□ 

□ 

Human Controlled 
(inploce 1948) 

O 

O 

O 

Humon Controlled 


• 

(introduced since 194 8) 

o 

o 

o 

Noturol/Scmi-Nat jral 

Figure 6 . LAMP 

COVER CONNECTIVITY GRAPHS 

(Hcccntly added land 

covers arc poorly connected to traditional cover types! 
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.TABLE I. AERIAL PHOTOGRAPHIC DATA ELEMENT FILE FORMAT 


Column 


Data Element 


1-2 

3-4 

5-7 

e-io 

11-12 

13-17 


18-19 

20-24 

25 

26-28 


29 

30-32 


33 

34-36 

37-39 

40-41 

42-43 

44-45 

46-47 

48-49 


50-51 
52-53 
54-55 
56-57 
58-59 
60-61 
62-63 
, 64-65 


blank 

year photo mission was flown 

flight line of frame on which grid square being interpreted appears 

frame number on which grid square being interpreted appears 

blank 

Universal Transverse Mercator (UTM)’ grid coordinate row number for 
the square being interpreted NOTE: This coordinate is always the 

coordinate of the upper left hand corner of the grid square, 
blank 

UTM coordinate column number for the square being interpreted 
blank 

three digit land cover for land cover occupying greatest area of 
grid square (see Table 2) For any grid square, there must be a 
land cover in this entry, 
blank 

three digit land cover code for land cover taking up second most 
extensive area 1IOTE : If the entire square is taken up by one land 

cover, this entry will be zero (0) . 
blank 

three digit land cover code for land cover taking up thi rd most 

extensive area in grid square (This may also be zero.) 

number of occupied houses identifiable in t'e square (dwelling units) 

number of schools identified in the square 

number of churches identified in the square 

number of other government buildings in the square 

number of buildings associated with industrial activity in the 

square (includes oil tanks and storage areas) 

type of roads/trails in square (Choose only the highest category.) 

0 = no road or trail 

1 <= all-weather road 

2 «= usually passable to wheeled vehicles 

3 trail not normally usable lay wheeled vehicles (Choose only the 
number of bridges in square 

number of boats in square 

number of piers qr boat landing points in the square 

number of concrete drying platforms or basketball courts in square 

number of fish ponds in the square 

number of fish traps in the square 

number of marketplaces in the square 

number of commercial buildings in the square, including sari-sari 
store front houses where identifiable 


TABLE II. LAND COVER CLASSIFICATION FOR MARINDUQUE ISLAND, PHILIPPINES 


Computer Photo 

Coding No. Symbol Land Cover 


101 n 

Ar 

102 

Ac 

104 

Aom 

105 

Aoc 

108 

Ai 

109 

Aq 

112 c 

Am 

201° 

Fr 

202 b 

Fb 

203" 

Fg 

301° 

Rh 

302 

Rm 

303. 

503 bfe 

Rc 

Es 

60S; 

Pm 

V°«® 

Ta 

9 or' 

Ns 

902 

Hr 

903 

Wn 

904 

Wr 

906 b 

Wn 


Land intensively used for rice 

Corn, beans, forage crops, upland croplands 

Mixed orchard (dooryard orchards) 

Coconut orchard 
Fallow lands 

Aquatic lands (commercial fish ponds, salt beds, moss 
production, oyster beds) 

Mangrove swamps (including nipa palm areas) 

Natural forest stands 

Forest brushlnnd (go-back, parang ) and grassland ( aogon ) 
Natural bamboo (including all riverine tree vegetation) 
High density greater than 10 houses per centimeter on 
1:15,000 air photo 

Medium density 5 to 10 houses per centimeter on 1:15,000 
air photo 

Road strip and cross road community 
Strip mining 

Military base and governmental centers 

Airport 

Snnd/beach 

Rock, exposed bare earth 

Fresh water lakes and ponds greater than 7.5 hectares 
Permanent streams and river (including the entire width of 
river courses) 

Open salt water 


^100 series, agriculture 
Added land cover category for Marinduque study 
*j200 series, forest land 
e 300 series, residential land cover 
£500 series, extractive industry land cover 
600 series, public and semi-public land cover 
^700 series, transportation land cover 
"900 series, nonproductive land cover 


TABLE III. CHANGES IN TiiE NUMBER OF LAND COVER CATEGORIES 
MARINDUQUE ISLAND, PHILIPPINES, 1948-1972 


Land Cover Categories 



Dominant 3 

Secondary* 3 

Tertiary 0 

Total 

194 8 

9 

13 

14 

15 

1967" 

10 

14 

15 

18 

19 72 

15 

22 

22 

23 


^Dominant land cover — areally most extensive in a minimum of one 25 ha cell 
Secondary land cover -*• areally second most extensive in a minimum of one 
c 25 ha cell 

Tertiary land cover — areally the third most extensive in a minimum of one 
^25 ha cell 

1967 data for only the north central-northeastern portion of the Island 


VO VO VO 


TABLE 'IV. DOKINANx .uAND COVER 3H CELL C25 HA), K APINDUfiUE ISLAND, PHILIPPINES, 1943-1972 


Ponds 



Porest 

Cocor.~ 

P.ock/ 

Bare 


No. 

% ! 

j No. 

% 

No. 

64 

1.5 

952 

23.3 

0 

57 

1.4 

1,323 

33.0 

11 

6 

.5 

235 ] 

21.1 

0 

12 

1.1 

314 

23.1 

7 

6 

.5 

362 

32.4 

3 


1654.1 
1631 4 . 2 
SS j 7. 9 
39 3.0 
37 7.3 


No. 

% 

No. 


403 

10.1 

1,772 

44.2 

321 

3.0 

1,445 

36.1 

46 

4.1 

535 

52.4 

43 

4.3 

431 

43.1 

44 

3.9 

432 

33.7 


% No. 


353 
367 
107 
4 I Of 

1,3 1102 I 9.1 


Cnx-Scuare Results 


Period 


Value 


Significant Difference 


43 versus 1972 
43 versus 1572 

1967 versus 


Entire Island 

North central-northeastern portion 
Includes dooryard orchards and sugar cane 
Includes barboo thickets and brushland 

x *a = .01 [9] = 21-656 

x *a * . 01 [8] = 20-090 
X _ _ .1 ri * ■» “ 23.209 


114.3 

57.4 

21.2 

16.0 























































TABLE V. CYCLIC AND TREND CHANGES IN TRIFORMED LAND COVER UNITS, NORTE CENT RAL-NORTEEA STERN 

NARINDUQUE ISLAND , PHILIPPINES, 1948-1972 


'riformed Land Cover Units pension and Contraction Linear Regression Results 



Primary 

Land 

Cover 

Secondary 

Land 

Cover 

■Coconut 

Paddy 

Coconut 

Paddy 

Coconut 

Cogon 

Coconut 

Cogon 

Coconut 

Cogon 

Coconut 

Paddy 

Cog on 

Dry Crop 

Cogon 

Coconut 

Cogon 

Coconut 

Cocon 

Forest 

Ccgcn 

Forest 

Cogon 

None 

Cogon 

Dry C-op 

Cogon 

Cry Crop 

Salt Water 

None 

Salt Water 

Mangrove 


Land 

Cover 



1967 


No. 

Cells 


Non 
ogo 
Ncn 
add 
Dry Crop 
Dry Crop 
Forest 
Paddy 
Dry Crop 
None 

Dry Crop 
None 
None 
Coconut 
None 
None 


18 1.6 14 

13 3.3 34 

42 3.8 25 


40 
57 
29 
39 
34 

6 . 8 I 55 




a Data fitted to linear regression model of fora V « ax + b. rover unit 

»*•»*•« ’'“‘•w*. «» ■ ?«£! t »67.19. 1, 


Dependent variable ui = nraer u- --------- ----- lq , 7 1972-24 

Independent variable (X) = time (year data collected) with 1948=>1, 1967-19, 1972 2 

a = slope of the trend line through X dimension 

b = coint of intercept of trend line or. Y axis . . , 

r = decree of fit of actual data to mathematically determined equation (correlation 

r 2 =» Correlation* coefficient squared (Pearson product-moment correlation coefficient) 


I 

*1 





















*X Mann ;| 

a Standard Deviation s 

ia Percent Standard Deviation ;? 



